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In these TASI summer school lectures I discuss three topics in neutrino astro- 
physics: the solar neutrino problem, stellar cooling by neutrino emission, and the 
role of neutrinos in the nucleosynthesis that occurs within core-collapse supernovae. 



1 Introduction 

Part of the interest in neutrino astrophysics has to do with the fascinating 
interplay between nuclear and particle physics issues — e.g., whether neutri- 
nos are massive and undergo flavor oscillations, whether they have detectable 
electromagnetic moments, etc. — and astrophysical phenomena, such as the 
clustering of matter on large scales, the mechanisms responsible for the synthe- 
sis of nuclei, and the evolution of stars. The three lectures here are intended 
to illustrate this interplay. The first lecture reviews the solar neutrino problem 
which, along with the atmospheric neutrino problem, provides perhaps our 
strongest direct evidence that new physics lurks beyond the standard model. 
The second has to with the implications of neutrino properties — e.g., whether 
neutrinos are Dirac or Majorana particles — for stellar cooling. The final lec- 
ture describes the nucleosynthesis we think accompanies a supernova explosion, 
and why that synthesis is a delicate probe of neutrino oscillations. 

These lectures as well of those from several other 1998 TASI speakers share 
a common subtheme: how the extraordinary technical revolution in astronomy 
and astrophysics has made the microphysics of the universe more relevant. It is 
the precise data coming from the new generation of great observatories — maps 
of the cosmic microwave spectrum, precision measurements of the products 
of big bang nucleosynthesis, measurements of the solar neutrino spectrum, 
Hubble Space Telescope (HST) abundance distributions from early, metal- 
poor stars, detection of gamma ray bursts from cosmological sources — that 
allow us to form the connections between observations and the underlying 
microphysics. This is the driving force that is making the field of nuclear and 
particle astrophysics of such interest to both senior physicists and new students 
entering the field. 
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2 Solar Neutrinos 

More than three decades ago Ray Davis, Jr. and his collaborators! constructed 
a 0.615 kiloton C2CI4 radiochemical solar neutrino detector in the Homestake 
Gold Mine, one mile beneath Lead, South Dakota. Within a few years it was 
apparent that the number of neutrinos detected was considerably below the 
predictions of the standard solar model, that is, the standard theory of main 
sequence stellar evolution. 

Today the results from the 37 C1 detector, which have become quite accu- 
rate due to 30 years of careful measurement, have beenpaugmented by results 
from four other experiments, the SAGEB and GALLEXel gallium experiments 
and the KamiokandaB and SuperKamiokandeB water Cerenkov detectors. It 
now appears that the combined results are very difficult to explain — some 
have argued impossible — by any plausible change in the standard solar model 
(SSM). Thus most believe that the answer to the solar neutrino problem is new 
particle physics, most likely some effect like solar neutrino oscillations associ- 
ated with massive neutrinos. With the recent news that SuperKamiokande sees 
direct evidence for oscillations in the azimuthal dependence of atmosphericQ 
neutrinos, it seems that we may be on the threshold of a major discovery. 

The purpose of this first (and longest) lecture is to summarize the solar 
neutrino problem and the arguments that it represents new particle physics. 

2.1 The Standard Solar ModeM 

Solar models trace the evolution of the sun over the past 4.6 billion years of 
main sequence burning, thereby predicting the present-day temperature and 
composition profiles of the solar core that govern neutrino production. Stan- 
dard solar models share four basic assumptions: 

* The sun evolves in hydrostatic equilibrium, maintaining a local balance 
between the gravitational force and the pressure gradient. To describe this 
condition in detail, one must specify the equation of state as a function of 
temperature, density, and composition. 

* Energy is transported by radiation and convection. While the solar envelope 
is convective, radiative transport dominates in the core region where ther- 
monuclear reactions take place. The opacity depends sensitively on the solar 
composition, particularly the abundances of heavier elements. 

* Thermonuclear reaction chains generate solar energy. The standard model 
predicts that over 98% of this energy is produced from the pp chain conversion 
of four protons into 4 He (see Fig. 1) 

4p He + 2e+ + 2v c (1) 
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with proton burning through the CNO cycle contributing the remaining 2%. 
The sun is a large but slow reactor: the core temperature, T c <~ 1.5 • 10 7 K, 
results in typical center-of-mass energies for reacting particles of ~ 10 keV, 
much less than the Coulomb barriers inhibiting charged particle nuclear re- 
actions. Thus reaction cross sections are small: in most cases, as laboratory 
measurements are only possible at higher energies, cross section data must be 
extrapolated to the solar energies of interest. 

* The model is constrained to produce today's solar radius, mass, and lumi- 
nosity. An important assumption of the standard model is that the sun was 
highly convective, and therefore uniform in composition, when it first entered 
the main sequence. It is furthermore assumed that the surface abundances of 
metals (nuclei with A > 5) were undisturbed by the subsequent evolution, and 
thus provide a record of the initial solar metallicity. The remaining param- 
eter is the initial 4 He/H ratio, which is adjusted until the model reproduces 
the present solar luminosity after 4.6 billion years of evolution. The resulting 
4 Hc/H mass fraction ratio is typically 0.27 ± 0.01, which can be compared 
to the big-bang value of 0.23 ± 0.01. Note that the sun was formed from 
previously processed material. 

The model that emerges is an evolving sun. As the core's chemical com- 
position changes, the opacity and core temperature rise, producing a 44% 
luminosity increase since the onset of the main sequence. The temperature 
rise governs the competition between the three cycles of the pp chain: the 
ppl cycle dominates below about 1.6 TO 7 K; the ppll cycle between (1.7-2.3) 
•10 7 K; and the ppIII above 2.4 -10 7 K. The central core temperature of today's 
SSM is about 1.55 -10 7 K. 

The competition between the cycles determines the pattern of neutrino 
fluxes. Thus one consequence of the thermal evolution of our sun is that the 
8 B neutrino flux, the most temperature-dependent component, proves to be 
of relatively recent origin: the predicted flux increases exponentially with a 
doubling period of about 0.9 billion years. 

A final aspect of SSM evolution is the formation of composition gradients 
on nuclear burning timescales. Clearly there is a gradual enrichment of the 
solar core in 4 He, the ashes of the pp chain. Another element, 3 He, is a sort of 
catalyst for the pp chain, being produced and then consumed, and thus even- 
tually reaching some equilibrium abundance. The timescale for equilibrium 
to be established as well as the eventually equilibrium abundance are both 
sharply decreasing functions of temperature, and thus increasing functions of 
the distance from the center of the core. Thus a steep 3 He density gradient is 
established over time. 

The SSM has had some notable successes. From helioseismology the sound 
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Figure 1: The solar pp chain. 

speed profile c(r) has been very accurately determined for the outer 90% of the 
sun, and is in excellent agreement with the SSM. Such studies verify important 
predictions of the SSM, such as the depth of the convective zone. However the 
SSM is not a complete model in that it does not explain all features of solar 
structure, such as the depletion of surface Li by two orders of magnitude. This 
is usually attributed to convective processes that operated at some epoch in 
our sun's history, dredging Li to a depth where burning takes place. 

The principal neutrino-producing reactions of the pp chain and CNO cycle 
are summarized in Table 1. The first six reactions produce f3 decay neutrino 
spectra having allowed shapes with endpoints given by E™ ax . Deviations from 
an allowed spectrum occur for 8 B neutrinos because the 8 Be final state is 
a broad resonance. The last two reactions produce line sources of electron 
capture neutrinos, with widths ~ 2 keV characteristic of the temperature of 
the solar core. Measurements of the pp, Be, and 8 B neutrino fluxes will 
determine the relative contributions of the ppl, ppll, and ppIII cycles to solar 
energy generation. As discussed above, and as later illustrations will show 
more clearly, this competition is governed in large classes of solar models by 
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Table 1: Solar neutrino sources and the flux predictions of the BP98 and Brun/Turck- 
Chieze/Morel SSMs in ern _2 s _1 . 



Source 


W nax ( MeV ) 
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BTCM98 


P + P 


-* 2 H + c+ 


+ V 


0.42 


5.94E10 


5.98E10 


13 N _ 


> 13 C + e+ - 


V v 
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> 15 N + c+ + v 
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18.77 
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7 Be + 


e" ^ 7 Li -4 


- V 


0.86 (90%) 


4.80E9 


4.70E9 








0.38 (10%) 






p + e~ 


+ p ^ 2 H 


+ v 


1.44 


1.39E8 


1.41E8 



a single parameter, the central temperature T c . The -flux predictions of the 
1998 calculations of Bahcall, Basu, and Pinsonneault H (BP98) and of Brun, 
Turck-Chieze and Morel □ are included in Table 1. 

2.2 Solar Neutrino Detection^^ 

Let us start with a brief reminder about low energy neutrino-nucleus interac- 
tions in detectors. Consider the charged current reaction 

v e + (A,Z)^e- + (A,Z+l) (2) 

Because the momentum transfer to the nucleus is very small for solar neutrinos, 
it can be neglected in the weak propagator, leading to an effective contact 
current-current interaction. If we begin with the simplest (though fictitious) 
case of the free neutron decay n — ► p, the corresponding transition amplitude 
is then 

G F 

Sfi = -^=cos6» c u(p)7m(1 ~ 9Al5)u{n)u(e)^{\ - ^)u{v) (3) 

where Gf is the weak coupling constant measured in muon decay and cos# c 
gives the amplitude for the weak interaction to connect the u quark to its first- 
generation partner, the d quark. The origin of this effective amplitude is the 
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underlying standard model predictions for the elementary quark and lepton 
currents. The weak interactions at this level are predicted by the standard 
model to be exactly left handed. Experiment shows that the effective coupling 
of the W boson to the nuclcon is governed by 7^(1 — 5.475), as noted above, 
where gA ~ 1.26. The axial coupling is thus shifted from its underlying value 
by the strong interactions responsible for the binding of the quarks within the 
nucleon. 

If an isolated nucleon were the target, one could proceed to calculate the 
cross section from the effective nucleon current given above. The extension 
to nuclear systems traditionally begins with the observation that nucleons in 
the nucleus are rather nonrelativistic, v/c ~ 0.1. The amplitude u(p)7 ,i (l — 
gA"i5)u(n) can be expanded in powers of p/M . The leading vector and axial 
operators are readily found to be 

7o : 1 

7 : p/M <~ v/c 
7o7 5 : a ■ p/M ~ v/c 
775 : <? 

Thus it is the time-like part of the vector current and the space-like part of 
the axial- vector current that survive in the nonrelativistic limit. 

(In a nucleus these currents must be corrected for the presence of meson 
exchange contributions. The corrections to the vector charge and axial three- 
current, which we just pointed out survive in the nonrelativistic limit, are of 
order (v/c) 2 <~ 1%. Thus the naive one-body currents are a very good ap- 
proximation to the nuclear currents. In contrast, exchange current corrections 
to the axial charge and vector three-current operators are of order v/c, and 
thus of relative order 1. This difficulty for the vector three-current can be 
largely circumvented, because current conservation as embodied in the gener- 
alized Siegert's theorem allows one to rewrite important parts of this operator 
in terms of the vector charge operator. In the long-wavelength limit appro- 
priate to solar neutrinos, all terms unconstrained by current conservation do 
not survive. In effect, one has replaced a current operator with large two-body 
corrections by a charge operator with only small corrections. In contrast, the 
axial charge operator is significantly altered by exchange currents even for 
long-wavelength processes like (5 decay. Typical axial-charge (3 decay rates are 
enhanced by ~ 2 because of exchange currents.) 

If such a nonrelativistic reduction is done for our single current one obtains 

Sfi „ cos^%(^(p)^(n)u(e)7 (l - ts)«M 
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-${p)g A e4>{ri) ■ u(e)7(l - i<s)u{v)) 



(4) 



where the 4> are now two-component Pauli spinors for the nucleons. The above 
result can be generalized to include v e reactions by introducing the isospin 
operators t± where r + | n) = | p) and r_| p) = | n), with all other matrix 
elements being zero. Thus 

<fi (p)acf>(n) -» 4>\N)aT±(t>{N). 

This result easily generalizes to nuclear decay. Given our comments about 
exchange currents, the first step is the replacement 

A 

i=l 

C7T± -> ^V(«)t ± («). 

i=l 

Plugging 5/i into the standard cross section formula (which involves an average 
over initial and sum over final nuclear spins of the square of the transition 
amplitude) then yields the allowed nuclear matrix element 

A A 

^-^-(k/iiE^wiwi'+^k/iiE^^wiwi 2 )- ( 5 ) 

* i=l i=l 

Our initial calculation for the nucleon treated that particle as structureless. 
Implicitly we assumed that the momentum transfer is much smaller than the 
inverse nucleon size. If we take 10 MeV as a typical solar neutrino momentum 
transfer, these quantities would be in the ratio 1:20. For a light nucleus, the 
corresponding result might be 1:10. This long-wavelength approximation in 
combination with the nonrelativistic approximation yields the allowed result, 
where only Fermi and Gamow- Teller operators survive. These are the spin- 
independent and spin-dependent operators appearing above. 

The Fermi operator is proportional to the isospin raising/lowering opera- 
tor: in the limit of good isopsin, which typically is good to 5% or better in the 
description of low-lying nuclear states, it can only connect states in the same 
isospin multiplet, that is, states with a common spin-spatial structure. If the 
initial state has isospin (Tj, Mti), this final state has (Tj, Mt% ± 1) for (v, e~) 
and (P, e + ) reactions, respectively, and is called the isospin analog state (IAS). 
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In the limit of good isospin the sum rule for this operator in then particularly 
simple 

A A 

E^VtK/iiE t +«IN>I 2 = ^ViK /as iiE t +«h*>i 2 - \ N ~ Z \- ( fi ) 

/ 1 i=l 1 i=l 

The excitation energy of the IAS relative to the parent ground state can be 
estimated accurately from the Coulomb energy difference 

1 72&Z 

ElAS ~ ( 1.12AV3 + .78 - L293)MeV ' (?) 

The angular distribution of the outgoing electron for a pure Fermi (N, Z)+u — > 
(iV — 1, Z + 1) + e~ transition is 1 + (i cos ve , and thus forward peaked. Here 
(3 is the electron velocity. 

The Gamow- Teller (GT) response is more complicated, as the operator 
can connect the ground state to many states in the final nucleus. In general 
we do not have a precise probe of the nuclear GT response apart from weak 
interactions themselves. However a good approximate probe is provided by 
forward-angle (p,n) scattering off nuclei, a technique that has been developed 
in particular by experimentalists at the Indiana University Cyclotron Facility. 
The (p,n) reaction transfers isospin and thus is superficially like (V, e~). At 
forward angles (p,n) reactions involve negligible three-momentum transfers to 
the nucleus. Thus the nucleus should not be radially excited. It thus seems 
quite plausible that forward-angle (p,n) reactions probe the isospin and spin of 
the nucleus, the macroscopic quantum numbers, and thus the Fermi and GT 
responses. For typical transitions, the correspondence between (p,n) and the 
weak GT operators is believed to be accurate to about 10%. Of course, in a 
specific transition, much larger discrepancies can arise. 

The (p,n) studies demonstrate that the GT strength tends to concentrate 
in a broad resonance centered at a position S = Eqt — Eias relative to the 
IAS given by 

S~ (7.0 -28.9 ^ A Z ) MeV. (8) 

Thus while the peak of the GT resonance is substantially above the IAS for 
N ~ Z nuclei, it drops with increasing neutron excess. Thus 5 ~ for Pb. A 
typical value for the full width at half maximum r is ~ 5 MeV. 

The approximate Ikeda sum rule constrains the difference in the (v, e~) 
and {v,e + ) strengths 

E(l^&(". Ol 2 - \Mgr(u, e + )| 2 ) = 3(N Z) (9) 
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where 

1 A 

\Mg T (v,e-)\ 2 = ^—^KfW^^r+im)] 2 . (10) 

i— 1 

In many cases of interest in heavy nuclei, the strength in the (P, e + ) direction 
is largely blocked. For example, in a naive 2s\d shell model description of 37 C1, 
the p — > n direction is blocked by the closed neutron shell at N=20. Thus this 
relation can provide an estimate of the total [3~ strength. Experiment shows 
that the (3~ strength found in and below the GT resonance does not saturate 
the Ikeda sum rule, typically accounting for ~ (60 — 70) % of the total. Mea- 
sured and shell model predictions of individual GT transition strengths tend 
to differ systematically by about the same factor. Presumably the missing 
strength is spread over a broad interval of energies above the GT resonance. 
This is not unexpected if one keeps in mind that the shell model is an approxi- 
mate effective theory designed to describe the long wavelength modes of nuclei: 
such a model should require effective operators, renormalizecLfxom their bare 
values. Phenomenologically, the shell model seems to required g e /^ ~ 1.0 as 
well as a small spin-tensor term {a <g> Y2(f))j—i of relative strength ~ 0.1. 

The angular distribution of GT (N, Z) + v e -> (N- 1, Z+l) + e" reactions 
is 3 — /3 cos 8 '„ e , corresponding to a gentle peaking in the backward direction. 

The above discussion of allowed responses can be repeated for neutral cur- 
rent processes such as (u, v'). The analog of the Fermi operator contributes 
only to elastic processes, where the standard model nuclear weak charge is 
approximately the neutron number. As this operator does not generate tran- 
sitions, it is not yet of much interest for solar or supernova neutrino detection, 
though there are efforts to develop low-threshold detectors (e.g., cryogenic 
technologies) for recording the modest nuclear recoil energies. The analog of 
the GT response involves 

|M^^^| 2 -^ll(/llE^)^IN)| 2 - (11) 

i— 1 

The operator appearing in this expression is familiar from magnetic moments 
and magnetic transitions, where the large isovector magnetic moment (fi v ~ 
4.706) often leads to it dominating the orbital and isoscalar spin operators. 

Finally, there is one purely leptonic reaction of great interest, since it is the 
reaction exploited by Kamiokande and SuperKamiokande. Electron neutrinos 
can scattered off electrons via both charged and neutral current reactions. 
The cross section calculation is straight forward and will not be repeated here. 
Two features of the result are of importance for our later discussions, however. 
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Because of the neutral current contribution, heavy-flavor [v^ and v T ) also 
scatter off electrons, but with a cross section reduced by about a factor of 
seven at low energies. Second, for neutrino energies well above the electron 
rest mass, the scattering is sharply forward peaked. Thus this reaction allows 
one to exploit the position of the sun in separating the solar neutrino signal 
from a large but isotropic background. 

As I mentioned earlier, the first experiment performed was one exploiting 
the reaction 

37 Cl(^,e-) 37 Ar. 

As the threshold for this reaction is 0.814 MeV, the important neutrino sources 
are the 7 Be and 8 B reactions. The 7 Be neutrinos excite just the GT transition 
to the ground state, the strength of which is known from the electron capture 
lifetime of 37 Ar. The 8 B neutrinos can excite all bound states in 37 Ar, including 
the dominant transition to the IAS residing at an excitation of 4.99 MeV. The 
strength of excite-state GT transitions can be determined from the (3 decay 
37 Ca(/3 + ) 37 K, which is the isospin mirror reaction to 37 C1(^, e~) 37 Ar. The 
net result is that, for SSM fluxes, 78% of the capture rate should be due 
to 8 B neutrinos, and 15% to 7 Be neutrinos. The measured capture rate EJ 
2.56 ±0.16 ± 0.16 SNU (1 SNU = 10~ 36 capture/atom/sec) is about 1/3 the 
standard model value. 

Similar radiochemical experiments were done by the SAGE and GALLEX 
collaborations using a different target, 71 Ga. The special properties of this tar- 
get include its low threshold and an unusually strong transition to the ground 
state of 71 Ge, leading to a large pp neutrino cross section (see Fig. 2). The 
experimental capture rates are 66 ± 13 ± 6 and 76 ± 8 SNU for the SAGE and 
GALLEX detectors, respectively. The SSM prediction is about 130 SNulii 
Most important, since the pp flux is directly constrained by the solar lumi- 
nosity in all steady-state models, there is a minimum theoretical value for the 
capture rate of 79 SNU, given standard model weak interaction physics. Note 
there are substantial uncertainties in the 71 Ga cross section due to 7 Be neutrino 
capture to two excited states of unknown strength. These uncertainties were 
greatly reduced by direct calibrations of both detectors using 51 Cr neutrino 
sources. 

The remaining experiments, Kamiokande II /III and SuperKamiokande, 
exploited water Cerenkov detectors to view solar neutrinos event-by-event. 
Solar neutrinos scatter off electrons, with the recoiling electrons producing the 
Cerenkov radiation that is then recorded in surrounding phototubes. Thresh- 
olds are determined by background rates; SuperKamiokande is currently op- 
erating with a trigger at approximately six MeV. The initial experiment, 
Kamiokande II/III, found a flux of 8 B neutrinos of (2.91 ±0.24 ± 0.35) • 
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Figure 2: Level scheme for 71 Ge showing the excited states that contribute to absorption of 
pp, 7 Be, 51 Crm and 8 B neutrinos. 



10 6 /cm 2 s after about a decade of measurement. Its much larger successor Su- 
pcrKamiokande, with a 22.5 kiloton fiducial volume, yielded the result (2.37 ± 
0.06 ± 0.08) • 10 6 /cm 2 s after the first 374 days of measurements. This is about 
36% of the SSM flux. This result continues to improve in accuracy. 



2.3 Uncertainties in Standard Solar Model Parameters 

The pattern of solar neutrino fluxes that has emerged from these experiments 
is 

0(pp)~ 0.9 4> ssm (pp) 
0( 7 Be) - 
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0( S B) ~ O.43 bbM ( 8 B). (12) 

A reduced 8 B neutrino flux can be produced by lowering the central temper- 
ature of the sun somewhat, as </>( 8 B)^ T c 18 . However, such an adjustment, 
either by varying the parameters of the SSM or by adopting some nonstan- 
dard physics, tends to push the 0( 7 Be)/0( 8 B) ratio to higher values rather 
than the low one of eq. (12), 



<K 7 Bc) 10 
0(8B) ~ c ■ 



(13) 



Thus the observations seem difficult to reconcile with plausible solar model 
variations: one observable (0( 8 B)) requires a cooler core while a second, the 
ratio </>( 7 Be)/0( 8 B), requires a hotter one. 

An initial question is whether this problem remains significant when one 
takes into account known uncertainties in the parameters of the SSM. While 
a detailed summary of standard model uncertainties would take us well be- 
yond the limits of these lectures, a qualitative discussion of pp chain nuclear 
uncertainties is appropriate. This nuclear microphysics has been the focus of 
a great deal of experimental work. The pp chain involves a series of non- 
resonant charged-particle reactions occurring at center-of-mass energies that 
are well below the height of the inhibiting Coulomb barriers. As the resulting 
small cross sections preclude laboratory measurements at the relevant energies, 
one must extrapolate higher energy measurements to threshold to obtain solar 
cross sections. This extrapolation is often discussed in terms of the astrophys- 
ical S-factor 

a(E) = ^-eM-^v) (14) 

where r\ = Zl ^ 2a ; with a the fine structure constant and (3 = v/c the relative 
velocity of the colliding particles. This parameterization removes the gross 
Coulomb effects associated with the s-wave interactions of charged, point-like 
particles. The remaining energy dependence of S(E) is gentle and can be ex- 
pressed as a low-order polynomial in E. Usually the variation of S(E) with E is 
taken from a direct reaction model and then used to extrapolate higher energy 
measurements to threshold. The model accounts for finite nuclear size effects, 
strong interaction effects, contributions from other partial waves, etc. As labo- 
ratory measurements are made with atomic nuclei while conditions in the solar 
core guarantee the complete ionization of light nuclei, additional corrections 
must be made to account for the different electronic screening environments. 

Recently a large working group met at a workshop sponsored by the Insti- 
tute for Nuclear Theory, University of Washington, to review past work on the 
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nuclear reactions of the pp chain and CNO cycle, to recommend best values 
and appropriate errors, and to identify specific issues in experiment and theory 
where additional work is needed. The results will soon be published in Reviews 
of Modern Physics. I will not attempt a summary here, but will give one or 
two highlights. 

The most significant recommend change involves_the reaction 7 Be(p, 7) 8 B, 
where the standard Si7(0)~ 22.4 eVb is that givenliJ by Johnson et al. Mea- 
surements of Si7(E) are complicated by the need to use radioactive targets and 
thus to determine the areal density of the 7 Be target nuclei. Two techniques 
have been employed, measuring the rate of 478 keV photons from 7 Be decay or 
counting the daughter 7 Li nuclei via the reaction 7 Li (d,p) 8 Li. The low-energy 
data sets for Si7(E) disagree by 25%. This is a systematic normalization prob- 
lem as each data set is consistent with theory in its dependence on E. The 
energy dependence below ~ 500 keV is believed to be quite simple as it is 
determined by the asymptotic nuclear wave function. 

The Seattle working groruxon Si7(E) found that only one low-energy data 
set, that of Filippone et al. ta, was described in the published literature in 
sufficient detail to be evaluated. The target activity in that experiment had 
been measured by both 478 keV gamma rays and by the (d,p) reaction, with 
consistent results. The resulting recommended value was thus based on this 
measurement, yielding 

5i 7 (0) = 19±^eV b, la. (15) 

The 3 He(a, 7) 7 Be reaction has been measured by two techniques, by count- 
ing the capture 7 rays and by detecting the resulting 7 Be activity. While the 
two techniques have been used by several groups and have yielded separately 
consistent results, the capture 7 ray value Si7(0) = 0.507 ± 0.016 keV b is not 
in good agreement with the 7 Be activity value 0.572 ± 0.026 keV-b. The Seat- 
tle working group concluded that the evidence for a systematic discrepancy of 
unknown origin was reasonably strong and recommended that standard proce- 
dures be used in assigning a suitably expanded error. The recommended value 
S 34 (0) is 0.53 ± 0.05. 

These and other recommended values were recently incorporated into the 
BP98 solar model calculation. While the workshop's recommended values in- 
volve no qualitative changes, there is some broadening of error bars. The 
downward shift in Si7(0) leads to a lower 8 B flux. The workshop's Reviews of 
Modern Physics article summarizes a substantial amount of work on topics not 
discussed here: screening effects, weak radiative corrections to and exchange 
current effects on p+p, the atomic physics of 7 Be + e~, etc. Much of this dis- 
cussion was useful in evaluating possible uncertainties in solar microphysics, 
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and in identifying opportunities for reducing these uncertainitics. 

Are uncertainties in the parameters of the SSM a significant source of 
uncertainty? The S-factors discussed above comprise one set of parameters, 
but there are others: the solar lifetime, the opacities, the solar luminosity, etc. 
In order to answer this question while also taking into account correlations 
among the fluxes when input parameters are varied, first Bahcall and Ulrichli3 
and later Bahcall and HaxtonEZI constructed 1000 SSMs by randomly varying 
five input parameters, the primordial heavy-element-to-hydrogen ratio Z/X 
and S(0) for the p-p, 3 He- 3 He, 3 He- 4 He, and p- 7 Be reactions, assuming for each 
parameter a normal distribution with the mean and standard deviation. (These 
were the parameters assigned the largest uncertainties.) Smaller uncertainties 
from radiative opacities, the solar luminosity, and the solar age were folded 
into the results of the model calculations perturbatively. 

The resulting pattern of 7 Be and 8 B flux predictions is shown in Fig. 3. 
The elongated error ellipses indicate that the fluxes are strongly correlated. 
Those variations producing </>( 8 B) below 0.8</> SSM ( 8 B) tend to produce a re- 
duced 0( 7 Be), but the reduction is always less than 0.8. Thus a greatly reduced 
</>( 7 Be) cannot be achieved within the uncertainties assigned to parameters in 
the SSM. 

A similar exploration, but including parameter variations very far from 
their preferred values, was carried out by Castellani et al. tLa, who displayed 
their results as a function of the resulting core temperature T c . The pattern 
that emerges is striking (see Fig. 4) : parameter variations producing the same 
value of T c produce remarkably similar fluxes. Thus T c provides an excellent 
one-parameter description of standard model perturbations. Figure 4 also 
illustrates the difficulty of producing a low ratio of </>( 7 Be)/0( 8 B) when T c is 
reduced. 

The 1000-solar-model variations were made under the constraint of repro- 
ducing the solar luminosity. Those variations show a similar strong correlation 
with T c 

4>{pp) oc T- 1 - 2 0( 7 Be) cx T 8 0( 8 B) cx T c 18 . (16) 

Figures 3 and 4 offer a strong argument that reasonable variations in the pa- 
rameters of the SSM, or nonstandard changes in quantities like the metallicity, 
opacities, or solar age, cannot produce the pattern of fluxes deduced from 
experiment (eq. (12)). This would seem to limit possible solutions to errors 
either in the underlying physics of the SSM or in our understanding of neutrino 
properties. 
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Figure 3: SSM 7 Be and S B flux predictions. The dots represent the results of SSM calcu- 
lations where the input parameters were varied according to their assigned uncertainties, as 
described in the text. The 90% and 99% confidence level error ellipses are shown. 



2.4 Nonstandard Solar Models 

Nonstandard solar models include both variations of SSM parameters far out- 
side the ranges that are generally believed to be reasonable (some examples 
of which are given in Figure 4), and changes in the underlying physics of the 
model. The solar neutrino problem has been a major stimulus to models: in 
fact, most suggestions were motivated by the hope of producing a cooler sun 
(T c ~ 0.95T C ) that would avoid conflict with the results of the 37 C1 experi- 
ment. The suggestions included models with low heavy element abundances 
("low Z" models), in which one abandons the SSM assumption that the initial 
heavy element abundances are those we measure today at the sun's surface; 
periodically mixed solar cores; models where hydrogen is continually mixed 
into the core by turbulent diffusion or by convective mixing; and models where 
the solar core is partially supported by a strong central magnetic field or by its 
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Figure 4: The responses of the pp, 7 Be, and 8 B neutrino fluxes to the indicated variations in 
solar model input parameters, displayed as a function of the resulting central temperature 
T c . From Castellani et al. 



rapid rotation, thereby relaxing the SSM assumption that hydrostatic equilib- 
rium is achieved only .through the gas pressure gradient. A larger list is given 
by Bahcall and Davis E3. To illustrate the kinds of consequences such models 
have, two of these suggestions are discussed in more detail below. 

In low-Z models one postulates a reduction in the core metallicity from Z 
~ 0.02 to Z ~ 0.002. This lowers the core opacity (primarily because metals 
are very important to free-bound electron transitions), thus reducing T c and 
weakening the ppll and ppIII cycles. The attractiveness of low-Z models is 
due in part to the existence of mechanisms for adding heavier elements to 
the sun's surface. These include the infall of comets and other debris, as 
well as the accumulation of dust as the sun passes through interstellar clouds. 
However, the increased radiative energy transport in low-Z models leads to a 
thin convective envelope, in contradiction to interpretations of the 5-minutc 
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solar surface oscillations. A low He mass fraction also results. As diffusion of 
material from a thin convective envelope into the interior would deplete heavy 
elements at the surface, investigators have also questioned whether present 
abundances could have accumulated in low-Z models. Finally, the general 
consistency of solar heavy element abundances with those observed in other 
main sequence stars makes the model appear contrived. 

Models in which the solar core (~ 0.2 M©) is intermittently mixed break 
the standard model assumption of a steady-state sun: for a period of sev- 
eral million years (the thermal relaxation time for the core) following mixing, 
the usual relationship between the observed surface luminosity and rate of 
energy (and neutrino) production is altered as the sun burns out of equilib- 
rium. Calculations show that both the luminosity and the 8 B neutrino flux 
are suppressed while the sun relaxes back to the steady state. Such models 
have been considered seriously because of instabilities associated with large 
gradients in the 3 He abundance, which in equilibrium varies as ~ T~ e , where 
T is the local temperature. The resulting steep profile is unstable under finite 
amplitude displacements of a volume to smaller r: the energy released by the 
increased 3 He burning at higher T can exceed the energy in the perturbation. 
For a discussion of the plausibility of such_a trigger for core mixing, one can 
see the original work of Dilke and Gough E2l as well as a more recent critique 
by MerryfieldEi The possibility that continuous mixing on time scales of 3 He 
mixing could produce a flux pattern close to that observed (e.g., a suppression 
in both the 8 B neutrino flu* and the 7 Be/ 8 B flux ratio) was recently discussed 
by Cumming and HaxtonE3. 

This discussion of two of the more seriously explored nonstandard model 
possibilities illustrates how changes motivated by the solar neutrino problem 
often produce other, unwanted consequences. In particular, many experts feel 
that the good SSM agreement with helioseismology is likely to be destroyed 
by changes such as those discussed above. 

Figure 5 is an illustration by Hata et al.S of the flux predictions of several 
nonstandard models, including a low-Z model consistent with the 37 C1 results. 
As in the Castellani et al. exploration, the results cluster along a track that 
defines the naive T c dependence of the 0( 7 Be)/</>( 8 B) ratio, well separated from 
the experimental contours. 

There is now a popular argument that no such nonstandard model can 
solve the solar neutrino problem: if one assumes undistorted neutrino spectra, 
no combination of pp, 7 Be, and 8 B neutrino fluxes fits the experimental results 
wellEj. In fact, in an unconstrained fit, the required 7 Be flux is unphysical, 
negative by almost 3er. Thus, barring some unfortunate experimental error, it 
appears we are forced to look elsewhere for a solution. 
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(K 7 Be)/(K 7 Be) SSM 

Figure 5: The fluxes allowed by the combined results of the various solar neutrino experi- 
ments compared to the results of SSM variations and various nonstandard solar models. The 
solid line in the naive T c power law discussed in the text. From Hata et al. 

If experimental error, SSM parameter uncertainties, and nonstandard solar 
physics are ruled out as potential solutions, new particle physics is left as the 
leading possibility. Suggested particle physics solutions of the solar neutrino 
problem include neutrino oscillations, neutrino decay, neutrino magnetic mo- 
ments, and weakly interacting massive particles. Among these, the Mikhcycv- 
Smirnov-Wolfcnstein effect — neutrino oscillations enhanced by matter inter- 
actions — is widely regarded as the most plausible. 

2.5 Neutrino Oscillations 

One odd feature of particle physics is that neutrinos, which are not required by 
any symmetry to be massless, nevertheless must be much lighter than any of 
the other known fermions. For instance, the current limit on the v e mass is < 5 
eV. The standard model requires neutrinos to be massless, but the reasons are 
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not fundamental. Dirac mass terms mo, analogous to the mass terms for other 
fermions, cannot be constructed because the model contains no right-handed 
neutrino fields. Neutrinos can also have Majorana mass terms 

vlm L v L and ~v^m R v R (17) 

where the subscripts L and R denote left- and right-handed projections of 
the neutrino field v, and the superscript c denotes charge conjugation. The 
first term above is constructed from left-handed fields, but can only arise as 
a nonrenormalizable effective interaction when one is constrained to generate 
mi with the doublet scalar field of the standard model. The second term is 
absent from the standard model because there are no right-handed neutrino 
fields. 

None of these standard model arguments carries over to the more general, 
unified theories that theorists believe will supplant the standard model. In 
the enlarged multiplets of extended models it is natural to characterize the 
fermions of a single family, e.g., v e , e, u, d, by the same mass scale mo- Small 
neutrino masses are then frequently explained as a result of the Majorana 
neutrino masses. In the seesaw mechanism, 

M v ~(° T mD ). (18) 

Diagonalization of the mass matrix produces one light neutrino, mught ~ 
and one unobservably heavy, rrihoavy ~ mil- The factor (mu/mij) is the needed 
small parameter that accounts for the distinct scale of neutrino masses. The 
masses for the v e ,^n, and v T are then related to the squares of the correspond- 
ing quark masses m u , m c , and m t . Taking tiir ~ 10 16 GeV, a typical grand 
unification scale for models built on groups like SO (10), the seesaw mechanism 
gives the crude relation 

m Ve : m Vlx : <-» 2 • 1CT 12 : 2 • 1(T 7 : 3 • l(T 3 cV. (19) 

The fact that solar neutrino experiments can probe small neutrino masses, and 
thus provide insight into possible new mass scales that are far beyond the 
reach of direct accelerator measurements, has been an important theme of the 
field. 

Now one of the most interesting possibilities for solving the solar neutrino 
problem has to do with neutrino masses. For simplicity we will discuss just two 
neutrinos. If a neutrino has a mass m, we mean that as it propagates through 
free space, its energy and momentum are related in the usual way for this mass. 
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Thus if we have two neutrinos, we can label those neutrinos according to the 
eigenstates of the free Hamiltonian, that is, as mass eigenstates. 

But neutrinos are produced by the weak interaction. In this case, we have 
another set of eigenstates, the flavor eigenstates. We can define a v e as the 
neutrino that accompanies the positron in (3 decay. Likewise we label by 
the neutrino produced in muon decay. 

Now the question: are the eigenstates of the free Hamiltonian and of the 
weak interaction Hamiltonian identical? Most likely the answer is no: we know 
this is the case with the quarks, since the different families (the analog of the 
mass eigenstates) do interact through the weak interaction. That is, the up 
quark decays not only to the down quark, but also occasionally to the strange 
quark. (This is why we had a cos 9 C in our weak interaction amplitude: the 
amplitude for u — > s is proportional to sin0 c .) Thus we suspect that the weak 
interaction and mass eigenstates, while spanning the same two-neutrino space, 
are not coincident: the mass eigenstates \v\) and \v2) (with masses mi and 
7712 ) are related to the weak interaction eigenstates by 

\u e ) = cos0„|z/i) + sin6»„|^ 2 ) 

K) = -sm0„K) + cos6 v \v 2 ) (20) 

where 9 V is the (vacuum) mixing angle. 

An immediate consequence is that a state produced as a \v e ) or a \ at 
some time t — for example, a neutrino produced in f3 decay — does not remain 
a pure flavor eigenstate as it propagates away from the source. This is because 
the different mass eigenstates comprising the neutrino will accumulate differ- 
ent phases as they propagate downstream, a phenomenon known as vacuum 
oscillations (vacuum because the experiment is done in free space). To see 
the effect, suppose we produce a neutrino in some f3 decay where we measure 
the momentum of the initial nucleus, final nucleus, and positron. Thus the 
outgoing neutrino is a momentum eigenstatec3. At time t—Q 

\v{t = 0)) = \v e ) = cos0„|i/i) +sin6 v \v 2 ). (21) 

Each eigenstate subsequently propagates with a phase 

gi(k-x—uii) gi(/c'of— m^+fe 2 4) (22) 

But if the neutrino mass is small compared to the neutrino momentum/energy, 
one can write 

v /mf + fc2~fc(l + ^). (23) 
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Thus we conclude 

\v{t)) = e i (^- s - kt ^( m i+ m l) t / 4k ) 

x[cos9 v \ Vl )e lSm2t ' Ak + sin6 v \v 2 )e- tSm2t ' Ak ]. (24) 

We see there is a common average phase (which has no physical consequence) 
as well as a beat phase that depends on 

(5m 2 = to 2 — m\. (25) 

Now it is a simple matter to calculate the probability that our neutrino state 
remains a \u e ) at time t 

P Ve (t) = \(VeHt))\ 2 

= 1 - sin 2 26 v sin 2 ( | -» 1 - ^ sin 2 26 v (26) 

\ 4k J 2 

where the limit on the right is appropriate for large t. Now E ~ k, where E 
is the neutrino energy, by our assumption that the neutrino masses are small 
compared to k. Thus we can reinsert the units above to write the probability 
in terms of the distance x of the neutrino from its source, 

P l ,(x) = l-sin 2 2^sin 2 (^^). (27) 

(When one properly describes the neutrino state as a wave packet, the large- 
distance behavior follows from the eventual separation of the mass eigenstates.) 
If the the oscillation length 

£ ° = w <»> 

is comparable to or shorter than one astronomical unit, a reduction in the solar 
v e flux would be expected in terrestrial neutrino oscillations. 

The suggestion that the solar neutrino problem could be explained by 
neutrino oscillations was first made by Pontecorvo in 1958, who pointed out 
the analogy with Kq <-» Kq oscillations. From the point of view of particle 
physics, the sun is a marvelous neutrino source. The neutrinos travel a long 
distance and have low energies (~ 1 MeV), implying a sensitivity to 

5m 2 >\Q- 12 eV 2 . (29) 

In the seesaw mechanism, Sm 2 ~ m 2 ., so neutrino masses as low as m-i ~ 
1CP 6 cV could be probed. In contrast, terrestrial oscillation experiments with 
accelerator or reactor neutrinos are typically limited to 5m 2 > 0.1 eV 2 . 
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From the expressions above one expects vacuum oscillations to affect all 
neutrino species equally, if the oscillation length is small compared to an as- 
tronomical unit. This is somewhat in conflict with the data, as we have argued 
that the 7 Be neutrino flux is quite suppressed. Furthermore, there is a weak 
theoretical prejudice that 8 V should be small, like the Cabibbo angle. The first 
objection, however, can be circumvented in the case of "just so" oscillations 
where the oscillation length is comparable to one astronomical unit. In this 
case the oscillation probability becomes sharply energy dependent, and one 
can choose Sm 2 to preferentially suppress one component (e.g., the monochro- 
matic 7 Be neutrinos). This scenario has been explored by several groups and 
remains an interesting possibility. However, the requirement of large mixing 
angles remains. 

Below we will see that stars allow us to "get around" this problem with 
small mixing angles. In preparation for this, we first present the results above 
in a slightly more general way. The analog of eq. (24) for an initial muon 
neutrino (\v(t = 0)) = is 

\v{t)) = fJik-S-kt-iml+mDt/ik) 

x[-sin#> 1 )e t5m2t/4 ' £ +cos6 v \u 2 )e- lSm2t/4k ] (30) 

Now if we compare eqs. (24) and (30) we see that they are special cases of a 
more general problem. Suppose we write our initial neutrino wave function as 

\v(t = 0)) = a e (t = 0)K) + a^(t = 0)K). (31) 

Then eqs. (24) and (30) tell us that the subsequent propagation is described 
by changes in a e (x) and a^(x) according to (this takes a bit of algebra) 

. d ( a e \ _ 1 / -Sm 2 cos 29 v Sm 2 sin 29 v \ ( a e \ ,„„>. 

l dx \a^)~ 4E \ Sm 2 sin 29 v Sm 2 cos 2B V ) ) ' ^ ' 

Note that the common phase has been ignored: it can be absorbed into the 
overall phase of the coeeficients a e and a^, and thus has no consequence. Also, 
we have equated x = t, that is, set c = 1. 



2. 6 The Mikheyev-Smirnov- Wolfenstein Mechanism 

The view of neutrino oscillations changed when Mikheyev and Smirnov @ 
showed in 1985 that the density dependence of the neutrino effective mass, 
a phenomenon first discussed by Wolfenstein in 1978, could greatly enhance 
oscillation probabilities: a v e is adiabatically transformed into a as it tra- 
verses a critical density within the sun. It became clear that the sun was not 
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only an excellent neutrino source, but also a natural regenerator for cleverly 
enhancing the effects of flavor mixing. 

While the original work of Mikheyev and Smirnov was numerical, their 
phenomenon was soon understood analytically as a level-crossing problem. If 
one writes the neutrino wave function in matter as in eq. (31), the evolution 
of a e (x) and a^x) is governed by 

.d_ ( a e \ _ J_ / 2EV2G F p{x) - 6m 2 cos 26 v Sm 2 sin 26 v \ ( a e ' 
l dx\a fJ/ )~4:E \5m 2 sin28 v - 2E^/2G F p{x) + Sm 2 cos 20 v ) 

(33) 

where Gf is the weak coupling constant and p(x) the solar electron density. 
If p(x) = 0, this is exactly our previous result and can be trivially integrated 
to give the vacuum oscillation solutions of Sec. 2.5. The new contribution 
to the diagonal elements, 2E\/2Gfp{x), represents the effective contribution 
to M 2 that arises from neutrino-electron scattering. The indices of refraction 
of electron and muon neutrinos differ because the former scatter by charged 
and neutral currents, while the latter have only neutral current interactions. 
The difference in the forward scattering amplitudes determines the density- 
dependent splitting of the diagonal elements of the new matter equation. 

It is helpful to rewrite this equation in a basis consisting of the light and 
heavy local mass eigenstates (i.e., the states that diagonalize the right-hand 
side of the equation) , 

\v L (x)) = cos9(x)\v e ) - sin0(x)\v ll ) 
\v H {x)) = sm6(x)\is e ) +cos0(z)K>- ( 34 ) 
The local mixing angle is defined by 

sin 29 v 



sm2e(x) 



^X 2 (x)+sin 2 26 v 



cos2fl(x) = X ^ = (35) 

^X 2 (x) +sin 2 28 v 

where X(x) = 2\/2Gfp{x)E / 'Sm 2 — cos29 v - Thus 6(x) ranges from 6 V to tt/2 
as the density p(x) goes from to oo. 
If we define 

\v(x)) = a H (x)\v H (x)) + a L {x)\v L (x)), (36) 
the neutrino propagation can be rewritten in terms of the local mass eigenstates 

i—( aH \ = ( A(a: ) ia{x)\fa H \ , 3? . 
dx\dL J \—ia(x) -X(x) ) \cll ) 
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with the splitting of the local mass eigenstates determined by 

Sr™ 2 
~2E 

and with mixing of these eigenstates governed by the density gradient 



Stti 2 I 

2A(a?) = — \jx\x) + sin 2 29 v (38) 



E \ V2G F £p(x)sm26 v 



a(x) = [-=-*) , f rv 2 (39) 

v ' \5m 2 J X 2 (x) +sin 2 29 v 

The results above are quite interesting: the local mass eigenstates diagonalizc 
the matrix if the density is constant. In such a limit, the problem is no more 
complicated than our original vacuum oscillation case, although our mixing 
angle is changed because of the matter effects. But if the density is not con- 
stant, the mass eigenstates in fact evolve as the density changes. This is the 
crux of the MSW effect. Note that the splitting achieves its minimum value, 
sin 26 v , at a critical density p c = p(x c ) 

2y/2EG F p c = 6m 2 cos 20 v (40) 

that defines the point where the diagonal elements of the original flavor matrix 
cross. 

Our local-mass-eigenstate form of the propagation equation can be trivially 
integrated if the splitting of the diagonal elements is large compared to the off- 
diagonal elements, 



j(x) 



\(x) 



sin 2 29 v 8m 2 1 [X(x) 2 + sin 2 2<9„] 3 / 2 ^ ^ 



cos29 v 2E |J_M£)| sin 3 26>. 

1 p c dx 1 



a(x) 

a condition that becomes particularly stringent near the crossing point, 



, . s\tl 2 29 v 8m 2 1 , , A . 

7c = ifrc) = —7^ ——7- » I- (42) 



cos 26 v 2E 



J_ dp(x) 
p c dx \x—x c 



The resulting adiabatic electron neutrino survival probability 0, valid when 
7c > 1, is 

P, a c diab = X - + i cos 20 v cos 20i (43) 

where 9i = 9{xi) is the local mixing angle at the density where the neutrino 
was produced. 

The physical picture behind this derivation is illustrated in Fig. 6. One 
makes the usual assumption that, in vacuum, the v e is almost identical to the 
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Figure 6: Schematic illustration of the MSW crossing. The dashed lines correspond to the 
electron-electron and muon-muon diagonal elements of the M 2 matrix in the flavor basis. 
Their intersection defines the level-crossing density p c . The solid lines are the trajectories 
of the light and heavy local mass eigenstates. If the electron neutrino is produced at high 
density and propagates adiabatically, it will follow the heavy-mass trajectory, emerging from 
the sun as a v^. 

light mass eigenstate, ^l(O), i.e., m\ < m 2 and cos9 v ~ 1. But as the density 
increases, the matter effects make the v e heavier than the v^, with v e — > vh{x) 
as p{x) becomes large. The special property of the sun is that it produces u e s 
at high density that then propagate to the vacuum where they are measured. 
The adiabatic approximation tells us that if initially v e ~ i/jj (x) , the neutrino 
will remain on the heavy mass trajectory provided the density changes slowly. 
That is, if the solar density gradient is sufficiently gentle, the neutrino will 
emerge from the sun as the heavy vacuum eigenstate, ~ v^. This guarantees 
nearly complete conversion of z^ e s into u^s, producing a flux that cannot be 
detected by the Homestake or SAGE/GALLEX detectors. 

But this does not explain the curious pattern of partial flux suppressions 
coming from the various solar neutrino experiments. The key to this is the 
behavior when 7 C < 1. Our expression for 7(2;) shows that the critical region 
for nonadiabatic behavior occurs in a narrow region (for small 6 V ) surrounding 
the crossing point, and that this behavior is controlled by the derivative of the 
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density. This suggests an analytic strategy for handling nonadiabatic crossings: 
one can replace the true solar density by a simpler (integrable!) two-parameter 
form that is constrained to reproduce the true density and its derivative at 
the crossing point x c . Two convenient choices are the linear (p(x) — a + 
bx) and exponential (p(x) = ae~ bx ) profiles. As the density derivative at x c 
governs the nonadiabatic behavior, this procedure should provide an accurate 
description of the hopping probability between the local mass eigenstates when 
the neutrino traverses the crossing point. The initial and ending points Xi and 
Xf for the artificial profile are then chosen so that p(xt) is the density where 
the neutrino was produced in the solar core and p(xj) = (the solar surface), 
as illustrated in in Fig. 7. Since the adiabatic result (p^diab-j Qi e p enc i s on iy on 
the local mixing angles at these points, this choice builds in that limit. But 
our original flavor-basis equation can then be integrated exactly for linear and 
exponential profiles, with the results given in terms of parabolic cylinder and 
Whittaker functions, respectively. 

That result can be simplified further by observing that the nonadiabatic 
region is generally confined to a narrow region around x c , away from the end- 
points Xi and Xf. We can then extend the artificial profile to x = ±00, as 
illustrated by the dashed lines in Fig. 7. As the neutrino propagates adiabati- 
cally in the unphysical region x < Xi, the exact soluation in the physical region 
can be recovered by choosing the initial boundary conditions 

a^(— 00 J = — a^{— 00) = cost^e J -°° 

Off (—00) = o e (— 00) = sin0je l ^-°° x< - x ^ dx , (44) 

That is, 00)) will then adiabatically evolve to \v(xi)) = \v e ) as x goes from 
— 00 to Xi. The unphysical region x > Xf can be handled similarly. 

With some algebra a simple generalization of the adiabatic result emerges 
that is valid for all Sm 2 /E and 8 V 

P v? = \ + \ cos 26 v cos 20i(l - 2P Uop ) (45) 

where Phop is the Landau-Zener probability of hopping from the heavy mass 
trajectory to the light trajectory cm. traversing the crossing point. For the 
linear approximation to the densityEij'Ej, 

^op = e"^ /2 . (46) 

As it must by our construction, P Ve reduces to p^ dlab for 7 C 3> 1. When the 
crossing becomes nonadiabatic (e.g., j c <C 1 ), the hopping probability goes to 
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nonadiabatic 




Figure 7: The top figure illustrates, for one choice of sin 2 29 and 5m 2 , that the region 
of nonadiabatic propagation (solid line) is usually confined to a narrow region around the 
crossing point r c . In the lower figure, the solid lines represent the solar density and a 
linear approximation to that density that has the correct initial and final values, as well 
as the correct density and density derivative at r c . Thus the linear profile is a very good 
approximation to the sun in the vicinity of the crossing point. The MSW equations can 
be solved analytically for this wedge. By extending the wedge to ±oo (dotted lines) and 
assuming adiabatic propagation in these regions of unphysical density, one obtains the simple 
Landau-Zcncr result discussed in the text. 



1, allowing the neutrino to exit the sun on the light mass trajectory as a i/ e , 
i.e., no conversion occurs. 

Thus there are two conditions for strong conversion of solar neutrinos: 
there must be a level crossing (that is, the solar core density must be sufficient 
to render v e ~ vh{x{) when it is first produced) and the crossing must be 
adiabatic. The first condition requires that Sm 2 /E not be too large, and the 
second -f c > 1. The combination of these two constraints, illustrated in Fig. 

8, defines a triangle of interesting parameters in the sin 29 v plane, 

as Mikheyev and Smirnov found by numerically integration. A remarkable 
feature of this triangle is that strong v e — > conversion can occur for very 
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small mixing angles (sin 2 29 ~ 10 3 ), unlike the vacuum case. 




10" 4 10" 2 1 

sin 26 v 



Figure 8: MSW conversion for a neutrino produced at the sun's center. The upper shaded 
region indices thoses 5m 2 /E where the vacuum mass splitting is too great to be overcome 
by the solar density. Thus no level crossing occurs. The lower shaded region defines the 
region where the level crossing is nonadiabatic (7 C less than unity). The unshaded region 
corresponds to adiabatic level crossings where strong u e — ► will occur. 

One can envision superimposing on Fig. 8 the spectrum of solar neutri- 
nos, plotted as a function of for some choice of dm 2 . Since Davis sees 
some solar neutrinos, the solutions must correspond to the boundaries of the 

X 2 

triangle in Fig. 8. The horizontal boundary indicates the maximum 22g_ for 
which the sun's central density is sufficient to cause a level crossing. If a spec- 
trum properly straddles this boundary, we obtain a result consistent with the 
Homestake experiment in which low energy neutrinos (large 1/E) lie above 
the level-crossing boundary (and thus remain f e 's), but the high-energy neu- 
trinos (small 1/E) fall within the unshaded region where strong conversion 
takes place. Thus such a solution would mimic nonstandard solar models in 
that only the 8 B neutrino flux would be strongly suppressed. The diagonal 
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boundary separates the adiabatic and nonadiabatic regions. If the spectrum 
straddles this boundary, we obtain a second solution in which low energy neu- 
trinos lie within the conversion region, but the high-energy neutrinos (small 
1/E) lie below the conversion region and are characterized by 7 <C 1 at the 
crossing density. (Of course, the boundary is not a sharp one, but is char- 
acterized by the Landau-Zener exponential). Such a nonadiabatic solution is 
quite distinctive since the flux of pp neutrinos, which is strongly constrained 
in the standard solar model and in any steady-state nonstandard model by 
the solar luminosity, would now be sharply reduced. Finally, one can imag- 
ine "hybrid" solutions where the spectrum straddles both the level-crossing 
(horizontal) boundary and the adiabaticity (diagonal) boundary for small 9, 
thereby reducing the 7 Be neutrino flux more than either the pp or 8 B fluxes. 

What are the results of a careful search for MSW solutions satisfying 
the Homestakc, Kamiokandc/SuperKamiokande, and SAGE/GALLEX con- 
straints? This has been done by several groups: recent results will be discussed 
by N. Hata in his lectures. One solution, corresponding to a region surrounding 
Sm 2 <~ 6 • 10 _6 eV 2 and sin 2 29 v ~ 6 • 10~ 3 , is the hybrid case described above. 
It is commonly called the small-angle solution. A second, large-angle solution 
exists, corresponding to Sm 2 ~ 10~ 5 eV 2 and sin 2 29 v <~ 0.6. These solutions 
can be distinguished by their characteristic distortions of the solar neutrino 
spectrum. The survival probabilities P^ SW (E) for the small- and large- angle 
parameters given above are shown as a function of E in Fig. 9. 

The MSW mechanism provides a natural explanation for the pattern of 
observed solar neutrino fluxes. While it requires profound new physics, both 
massive neutrinos and neutrino mixing are expected in extended models. The 
small-angle solution corresponds to 6m 2 <~ 10~ 5 cV 2 , and thus is consistent 
with 1112 ~ few -10~ 3 eV. This is a typical v T mass in models where wir ~ 
togut- This mass is also reasonably close to atmospheric neutrino values. 
On the other hand, if it is the participating in the oscillation, this gives 
ma <~ 10 12 GeV and predicts a heavy v T ~ 10 eV. Such a mass is of great 
interest cosmologically as it would have consequences for supernova physics, 
the dark matter problem, and the formation of large-scale structure. 

2. 7 Spin-Flavor Oscillations 

If the MSW mechanism proves not to be the solution of the solar neutrino prob- 
lem, it still will have greatly enhanced the importance of solar neutrino physics: 
the existing experiments have ruled out large regions in the 5m 2 — sin 2 29 v plane 
(corresponding to nearly complete v e — > conversion) that remain hopelessly 
beyond the reach of accelerator neutrino oscillation experiments. 
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Figure 9: MSW survival probabilities P(E„) for typical small angle and large angle solutions. 



A number of other particle physics solutions have been considered, such as 
neutrino decay, matter-catalyzed neutrino decay, and solar energy transport 
by weakly interacting massive particles. But perhaps the most interesting 
possibility, apart from the MSW mechanism, was stimulated by suggestions 
that the 37 CI signal might be varying with a period comparable to the 11- 
year solar cycle. While the evidence for this has weakened, the original claims 
generated renewed interest in neutrino magnetic moment interactions with the 
solar magnetic field. 

The original suggestions by Cisneros and by Okun, Voloshyn, and Vysotsky 
envisioned the rotation 

Ve L -> v eR (47) 

producing a right-handed neutrino with sterile interactions in the standard 
model. With the discovery of the MSW mechanism, it was realized that matter 
effects would break the vacuum degeneracy of the v eL and v eR , suppressing the 
spin precession. Lim and MarcianoE3 and Akmedovcil pointed out that this 
difficulty was naturally circumvented for 
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as the different matter interactions of the v e and can compensate for the 
vacuum v e — mass difference, producing a crossing similar to the usual MSW 
mechanism. Such spin-flavor precession can then occur at full strength due to 
an off-diagonal (in flavor) magnetic moment. 

Quite relevant to this suggestion is the very strong limit on both diagonal 
and off-diagonal magnetic moments is imposed by studies of the red giant 
cooling process of plasmon decay into neutrinos 

7* -» v^j. (49) 

The result is < 3 • 10~ 12 /is, where /is is an electron Bohr magnetonEl. 
(This can be compared to a simple one- loop estimate il of the neutrino mag- 
netic moment of ~ 10~ 1s /ib, taking a typical dark matter value for the neutrino 
mass of a few eV.) If a magnetic moment at the red giant limit is assumed, it 
follows that solar magnetic field strengths of Bq > 10 6 G are needed to produce 
interesting effects. Since the location of the spin-flavor level crossing depends 
on the neutrino energy, such fields have to be extensive to affect an appreciable 
fraction of the neutrino spectrum. It is unclear whether these conditions can 
occur in the sun. This constraint from stellar cooling leads us naturally into 
our next lecture. 



3 Dirac and Majorana Neutrinos and Stellar Cooling 

3.1 The Neutrino Mass Matrix 

Consider a general 4n x 4n neutrino mass matrix where n is the number of 
flavors 
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\M* D Mr 

where each entry in this matrix is understood to be a n x n matrix operating in 
flavor space. The entries Mr, are the Dirac mass terms, while the Ml and Mr 
are the Majorana terms. The latter break the invariance of the Dirac equation 
under the transformation ip{x) — > e ia ip(x) associated with a conserved lepton 
number. Thus it is these terms that govern lepton-number-violating processes 
like double beta decay. 

One can proceed to diagonalize this matrix 
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U^Vi(x) with masses mj. 
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The eigenstates are two-component Majorana neutrinosu, yielding the proper 
2 x 2n = An degrees of freedom, where n is the number of flavors. We can 
recover the Majorana and Dirac limits: 

• If Mr = Ml = 0, the eigenstates of this matrix become pairwise degener- 
ate, allowing the 2n two-component eigenstates to be paired to form n four- 
component Dirac eigenstates. 

• If Mo — 0, the left- and right-handed components decouple, yielding n left- 
handed Majorana eigenstates with standard model interactions. 

There are interesting physical effects associated with these limits. Dirac 
neutrinos can have magnetic dipole, electric dipole (CP and T violating), and 
anapole (P violating) moments, as well as nonzero charge radii. Majorana 
neutrinos can have anapole moments but only transition magnetic and electric 
dipole moments. I mentioned in the previous lecture that transition magnetic 
moments were quite interesting in the context of MSW effects, as well as the 
use of Mr in the seesaw mechanism. 

3.2 Red Giants and Helium Burning 

We now consider the evolution off the main sequence of a solar-like star, with a 
mass above half a solar mass. As the hydrogen burning in the core progresses 
to the point that no more hydrogen is available, the stellar core consists of 
the ashes from this burning, 4 He. The star then goes through an interesting 
evolution: 

• With no further means of producing energy, the core slowly contracts, thereby 
increasing in temperature as gravity does work on the core. 

• Matter outside the core is still hydrogen rich, and can generate energy 
through hydrogen burning. Thus burning in this shell of material supports 
the outside layers of the star. Note as the core contracts, this matter outside 
the core also is pulled deeper into the gravitational potential. Furthermore, 
the shell H burning continually adds more mass to the core. This means the 
burning in the shell must intensify to generate the additional gas pressure to 
fight gravity. The shell also thickens as this happens, since more hydrogen is 
above the burning temperature. 

• The resulting increasing gas pressure causes the outer envelope of the star 
to expand by a larger factor, up to a factor of 50. The increase in radius 
more than compensates for the increased internal energy generation, so that a 
cooler surface results. The star reddens. Thus this class of star is named a red 
supergiant. 

• This evolution is relatively rapid, perhaps a few hundred million years: the 
dense core requires large energy production. The helium core is supported by 
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its degeneracy pressure, and is characterized by densities <~ 10 6 g/cm 3 . This 
stage ends when the core reaches densities and temperatures that allow helium 
burning through the reaction 

a + a + a^ 12 C + 7. (52) 

As this reaction is very temperature dependent (see below), the conditions for 
ignition are very sharply defined. This has the consequence that the core mass 
at the helium flash point is well determined. 

• The onset of helium burning produces a new source of support for the core. 
The energy release elevates the temperature and the core expands: He burn- 
ing, not electron degeneracy, now supports the core. The burning shell and 
envelope have moved outward, higher in the gravitational potential. Thus shell 
hydrogen burning slows (the shell cools) because less gas pressure is needed 
to satisfy hydrostatic equilibrium. All of this means the evolution of the star 
has now slowed: the red giant moves along the "horizontal branch" , as interior 
temperatures slowly elevate much as in the main sequence. 

The 3a process depends on some rather interesting nuclear physics. The 
first interesting "accident" involves the near degeneracy of the 8 Be ground 
state and two separated as: The 8 Be + ground state is just 92 keV above the 
2a threshold. The measured width of the 8 Be ground state is 2.5 eV, which 
corresponds to a lifetime of 

r m - 2.6 • l(T 16 s. (53) 

One can compare this number to the typical time for one a to pass another. 
The red giant core temperature is T 7 ~ 10 -> E ~ 8.6 keV. Thus v/c ~ 0.002. 
So the transit time is 



d 5/ 1 10~ 13 c 



cm 



v 0.002 3 • 10 10 cm/sec f 



8 • lQ- 2i s. (54) 



This is more than five orders of magnitude shorter than r m above. Thus when 
a 8 Be nucleus is produced, it lives for a substantial time compared to this naive 
estimate. 

To quantify this, we calculate the flux- averaged cross section assuming 
resonant capture 

(av) = (^f/^e- E ^T (55) 

where T is the 2a width of the 8 Be ground state. This is the cross section for 
the a + a reaction to form the compound nucleus then decay by a + a. But 
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since there is only one channel, this is clearly also the result for producing the 
compound nucleus 8 Bc. 

By multiplying the rate/volume for producing 8 Be by the lifetime of 8 Be, 
one gets the number of 8 Be nuclei per unit volume 

ATm s N a N a 
N(Be) = (av)T m 
1 + o aa 

N N 1 
1 + o aa 1 

= %^f?' 2 *- Er/kT - (56) 

Notice that the concentration is independent of Y. So a small T is not the 
reason we obtain a substantial buildup of 8 Be. This is easily seen: if the width 
is small, then the production rate of 8 Be goes down, but the lifetime of the 
nucleus once it is produced is longer. The two effects cancel to give the same 
8 Be concentration. One sees that the significant 8 Be concentration results from 
two effects: 1) a + a is the only open channel and 2) the resonance energy is 
low enough that some small fraction of the a + a reactions have the requisite 
energy. As E r = 92 keV, E r /kT = 10.67/T 8 (where T 8 is the temperature in 
1 8 K) so that 

N(Be) = Ar2 T -3/2 e -10.67/T 8 (Q 94 . (5?) 

So plugging in typical values of N a <~ 1.5T0 28 /cm 3 (corresponding to p a <~ 10 5 
g/cm 3 ) and T 8 ~ 1 yields 

^^^3.2x10 . (58) 

Salpeter suggested that this concentration would then allow a+ 8 Be^ 12 C 
to take place. Hoyle then argued that this reaction would not be fast enough 
to produce significant burning unless it was also resonant. Now the mass of 
8 Be + a, relative to 12 C,is 7.366 MeV, and each nucleus has ,V = 0+. Thus 
s-wave capture would require a + resonance in 12 C at ~ 7.4 MeV. No such 
state was then known, but a search by Cook, Fowler, Lauritsen, and Lauritsen 
revealed a + level at 7.644 MeV, with decay channels 8 Be+a and 7 decay to 
the 2+ 4.433 level in 12 C. The parameters are 

T Q - 8.9eV T 7 - 3.6 • l(T 3 cV. (59) 

The resonant cross section formula gives 

r 48 = AW^) 3/2 ^e-^ /fcT (60) 
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Plugging in our previous expression for iV( 8 Bc) yields 

r 48 = 7V3T 8 - 3 e- 42 - 9 / Ts (6.3 • l(r 54 cm 6 /s). (61) 
If we denote by cj 3q the decay rate of an a in our plasma, then 

u; 3a = 3A^ r -3 e -42.9/T 8(6 3 . 10 -54 cm 6 /sec) 

= ( 1.5^/c m 3 ^ 4 - 3 • 10 3 /sec)T 8 -3 e ---^. (62) 

Now the energy release per reaction is 7.27 MeV. Thus we can calculate 
the energy produced per gram, e: 

7.27McV 1.5 • 10 23 

e = w 3a 5 

3 g 

= (2.5 • ICPerg/g **)( - ■ J'T^-^/t. . (63) 

We can evaluate this at a temperature of T 8 <~ 1 to find 

e ~ (584ergs/g sec)( - -^j^ 3 ?■ ( 64 ) 

Typical values found in stellar calculations are in good agreement with this: 
typical red giant energy production is <~ 100 ergs per gram per second. 

To get a feel for the temperature sensitivity of this process, we can do a 
Taylor series expansion, finding 

<T) ~ (^) 40 Nl (65) 

J o 

This steep temperature dependence is the reason the He flash is delicately de- 
pendent on conditions in the core. 



3.3 Neutrino Magnetic Moments and He Ignition 

Prior to the helium flash, the degenerate He core radiates energy largely by 
neutrino pair emission. The process is the decay of a plasmon — which one 
can think of as a photon "dressed" by electron-hole excitations — thereby 
acquiring an effective mass of about 10 kcV. The photon couples to a neutrino 
pair through a electron particle-hole pair that then decays into aZ -t vv. 
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If this cooling is somehow enhanced, the degenerate helium core would 
be kept cooler, and would not ignite at the normal time. Instead it would 
continue to grow until it overcame the enhanced cooling to reach, once again, 
the ignition temperature. 

One possible mechanism for enhanced cooling is a neutrino magnetic mo- 
ment. Then the plasmon could directly couple to a neutrino pair. The strength 
of this coupling would depend on the size of the magnetic moment. 

A delay in the time of He ignition has several observable consequences, 
including changing the ratio of red giant to horizontal branch stars. Thus, using 
the standard theory of red giant evolution, investigators have attempted to 
determine what size of magnetic momcpt would produce unacceptable changes 
in the astronomy. The result is a limit E3 on the neutrino magnetic moment of 

fiij < 3 • 10~ 12 electron Bohr magnetons (66) 

as was mentioned earlier. This limit is more than two orders of magnitude 
more stringent than that from direct laboratory tests. 

This example is just one of a number of such constraints that can be 
extracted from similar stellar cooling arguments. The arguments above, for 
example, can be repeated for neutrino electric dipole moments. More inter- 
esting, it can be repeated for axion emission from red giants. Axions, the 
pseudoGolstone bosons associate with the solution of the strong CP problem 
suggested by Peccei and Quinn, are very light and can be produced radiatively 
within the red giant by the Compton process, by the Primakoff process off 
nuclei, or by emission from low-lying nuclear levels, such at from the 14 keV 
transition in 57 Fe. The net result is that axions of mass above a few eV are 
excluded; if axions have a direct coupling to electrons, so that the Compton 
process off electrons operates, the constraint is considerably tighter. 

A similar argument can be formulated for supernova cooling. During 
SN1987A the neutrino burst detected by 1MB and by Kamiokande was consis- 
tent with cooling on a timescale of about 4 seconds. Thus any process cooling 
the star more efficiently than neutrino emission would have shortened this time, 
while also reducing the flux in neutrinos. Large Dirac neutrino masses allow 
trapped neutrinos to scatter into sterile right-handed states. Right-handed 
neutrinos, lacking standard model interactions, would then escape the star 
(provided they do not scatter back into interacting left-handed states). Un- 
fortunately the upper bounds imposed on the neutrino mass are quite model 
dependent, ranging over (1-25) keV. 

The supernova cooling argument can also be repeated for axions. The 
window of sensitive runs from 1 eV (above this mass they are more strongly 
coupled than neutrinos, and thus cannot compete with neutrino cooling) to 
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about 0.01 eV (below this mass they are too weakly coupled to be produced 
on the timescale of supernova cooling). It is interesting that the supernova 
and red giant cooling limit on axions nearly meet: a small window may still 
exist around a few eV if the axion has no coupling to electrons. 

4 Supernovae, Supernova Neutrinos, and Nucleosynthesis 

Consider a massive star, in excess of 10 solar masses, burning the hydrogen in 
its core under the conditions of hydrostatic equilibrium. When the hydrogen 
is exhausted, the core contracts until the density and temperature are reached 
where 3a -^ 12 C can take place. The He is then burned to exhaustion. This 
pattern (fuel exhaustion, contraction, and ignition of the ashes of the previous 
burning cycle) repeats several times, leading finally to the explosive burning 
of 28 Si to Fe. For a heavy star, the evolution is rapid: the star has to work 
harder to maintain itself against its own gravity, and therefore consumes its 
fuel faster. A 25 solar mass star would go through all of these cycles in about 
7 My, with the final explosive Si burning stage taking a few days. The result 
is an "onion skin" structure of the precollapse star in which the star's history 
can be read by looking at the surface inward: there are concentric shells of H, 
4 He, 12 C, 16 and 20 Ne, 28 Si, and 56 Fe at the center. 

4-1 The Explosion Mechanism^ 

The source of energy for this evolution is nuclear binding energy. A plot of the 
nuclear binding energy S as a function of nuclear mass shows that the minimum 
is achieved at Fe. In a scale where the 12 C mass is picked as zero: 



12 c 


(5/nucleon = 


0.000 MeV 


16Q 


(5/nucleon = 


-0.296 MeV 


28 Si 


(5/nucleon = 


-0.768 MeV 


40 Ca 


(5/nucleon = 


-0.871 MeV 


56 Fe 


(5/nucleon = 


-1.082 MeV 


72 Ge 


(5/nucleon = 


-1.008 MeV 


98 Mo 


(5/nucleon = 


-0.899 Mev 



Thus once the Si burns to produce Fe, there is no further source of nuclear 
energy adequate to support the star. So as the last remnants of nuclear burning 
take place, the core is largely supported by degeneracy pressure, with the 
energy generation rate in the core being less than the stellar luminosity. The 
core density is about 2 xlO 9 g/cc and the temperature is kT ~ 0.5 MeV. 

Thus the collapse that begins with the end of Si burning is not halted 
by a new burning stage, but continues. As gravity does work on the matter, 
the collapse leads to a rapid heating and compression of the matter. As the 
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nucleons in Fe are bound by about 8 MeV, sufficient heating can release as and 
a few nucleons. At the same time, the electron chemical potential is increasing. 
This makes electron capture on nuclei and any free protons favorable, 



e +p 



v e + n. 



(67) 



Note that the chemical equilibrium condition is 



fl e + fi p = n n + (E„). 



(68) 



Thus the fact that neutrinos are not trapped plus the rise in the electron Fermi 
surface as the density increases, lead to increased neutronization of the matter. 
The escaping neutrinos carry off energy and lcpton number. Both the electron 
capture and the nuclear excitation and disassociation take energy out of the 
electron gas, which is the star's only source of support. This means that the 
collapse is very rapid. Numerical simulations find that the iron core of the star 
(<~ 1.2-1.5 solar mases) collapses at about 0.6 of the free fall velocity. 

In the early stages of the infall the v e s readily escape. But neutrinos 
are trapped when a density of ~ 10 12 g/cm 3 is reached. At this point the 
neutrinos begin to scatter off the matter through both charged current and 
coherent neutral current processes. The neutral current neutrino scattering 
off nuclei is particularly important, as the scattering cross section is off the 
total nuclear weak charge, which is approximately the neutron number. This 
process transfers very little energy because the mass energy of the nucleus is 
so much greater than the typical energy of the neutrinos. But momentum is 
exchanged. Thus the neutrino "random walks" out of the star. When the 
neutrino mean free path becomes sufficiently short, the "trapping time" of the 
neutrino begins to exceed the time scale for the collapse to be completed. This 
occurs at a density of about 10 12 g/cm 3 , or somewhat less than 1% of nuclear 
density. After this point, the energy released by further gravitational collapse 
and the star's remaining lcpton number are trapped within the star. 

If we take a neutron star of 1.4 solar masses and a radius of 10 km, an 
estimate of its binding energy is 



Thus this is roughly the trapped energy that will later be radiated in neutrinos. 

The trapped lepton fraction Yl is a crucial parameter in the explosion 
physics: a higher trapped Yl leads to a larger homologous core, a stronger 
shock wave, and easier passage of the shock wave through the outer core, as 
will be discussed below. Most of the lepton number loss of an infalling mass 
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element occurs as it passes through a narrow range of densities just before 
trapping. The reasons for this arc relatively simple: on dimensional grounds 
weak rates in a plasma go as T 5 , where T is the temperature. Thus the electron 
capture rapidly turns on as matter falls toward the trapping radius, and lepton 
number loss is maximal just prior to trapping. Inelastic neutrino reactions have 
an important effect on these losses, as the coherent trapping cross section goes 
as E 2 and is thus least effective for the lowest energy neutrinos. As these 
neutrinos escape, inelastic reactions repopulate the low energy states, allowing 
the neutrino emission to continue. 

The velocity of sound in matter rises with increasing density. The inner 
homologous core, with a mass M^c ~ 0-6 — 0.9 solar masses, is that part of 
the iron core where the sound velocity exceeds the infall velocity. This allows 
any pressure variations that may develop in the homologous core during infall 
to even out before the collapse is completed. As a result, the homologous core 
collapses as a unit, retaining its density profile. That is, if nothing were to 
happen to prevent it, the homologous core would collapse to a point. 

The collapse of the homologous core continues until nuclear densities are 
reached. As nuclear matter is rather incompressible (<~ 200 MeV/f 3 ), the 
nuclear equation of state is effective in halting the collapse: maximum densities 
of 3-4 times nuclear are reached, e.g., perhaps 6 • 10 14 g/cm 3 . The innermost 
shell of matter reaches this supernuclear density first, rebounds, sending a 
pressure wave out through the homologous core. This wave travels faster than 
the infalling matter, as the homologous core is characterized by a sound speed 
in excess of the infall speed. Subsequent shells follow. The resulting series of 
pressure waves collect near the sonic point (the edge of the homologous core). 
As this point reaches nuclear density and comes to rest, a shock wave breaks 
out and begins its traversal of the outer core. 

Initially the shock wave may carry an order of magnitude more energy 
than is needed to eject the mantle of the star (less than 10 51 ergs). But as 
the shock wave travels through the outer iron core, it heats and melts the iron 
that crosses the shock front, at a loss of ~ 8 MeV/nucleon. The enhanced 
electron capture that occurs off the free protons left in the wake of the shock, 
coupled with the sudden reduction of the neutrino opacity of the matter (recall 
o 'coherent ~ N 2 ), greatly accelerates neutrino emission. This is another energy 
loss. [Many numerical models predict a strong "breakout" burst of v e s in the 
few milliseconds required for the shock wave to travel from the edge of the 
homologous core to the neutrinosphere at p ~ 10 12 g/cm 3 and r ~ 50 km. 
The neutrinosphere is the term from the neutrino trapping radius, or surface 
of last scattering.] The summed losses from shock wave heating and neutrino 
emission are comparable to the initial energy carried by the shock wave. Thus 
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most numerical models fail to produce a successful "prompt" hydrodynamic 
explosion. 

Most of the attention in the past decade focused on two explosion scenarios. 
In the prompt mechanism described above, the shock wave is sufficiently strong 
to survive the passage of the outer iron core with enough energy to blow off 
the mantle of the star. The most favorable results were achieved with smaller 
stars (less than 15 solar masses) where there is less overlying iron, and with soft 
equations of state, which produce a more compact neutron star and thus lead 
to more energy release. In part because of the lepton number loss problems 
discussed earlier, now it is widely believed that this mechanism fails for all but 
unrealistically soft nuclear equations of state. 

The delayed mechanism begins with a failed hydrodynamic explosion; after 
about 0.01 seconds the shock wave stalls at a radius of 200-300 km. It exists in 
a sort of equilibrium, gaining energy from matter falling across the shock front, 
but loosing energy to the heating of that material. However, after perhaps 0.5 
seconds, the shock wave is revived due to neutrino heating of the nucleon 
"soup" left in the wake of the shock. This heating comes primarily from 
charged current reactions off the nucleons in that nucleon gas; quasielastic 
scattering also contributes. This high entropy radiation-dominated gas may 
reach two McV in temperature. The pressure exerted by this gas helps to push 
the shock outward. It is important to note that there are limits to how effective 
this neutrino energy transfer can be: if matter is too far from the core, the 
coupling to neutrinos is too weak to deposite significant energy If too close, 
the matter may be at a temperature (or soon reach a temperature) where 
neutrino emission cools the matter as fast or faster than neutrino absorption 
heats it. The term "gain radius" is used to describe the region where useful 
heating is done. 

This subject is still controversial and unclear. The problem is numerically 
challenging, forcing modelers to handle the difficult hydrodynamics of a shock 
wave; the complications of the nuclear equation of state at densities not yet 
accessible to experiment; modeling in two or three dimensions; handling the 
slow diffusion of neutrinos; etc. Not all of these aspects can be handled rea- 
sonably at the same time, even with existing supercomputers. Thus there is 
considerable disagreement about whether we have any supernova model that 
succeeds in ejecting the mantle. 

However the explosion proceeds, there is agreement that 99% of the 3 -10 53 
ergs released in the collapse is radiated in neutrinos of all flavors. The time 
scale over which the trapped neutrinos leak out of the protoneutron star is 
about three seconds. Through most of their migration out of the protoneutron 
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star, the neutrinos are in flavor equilibrium 

e.g., v e + v e f-> + v^. (70) 

As a result, there is an approximate equipartition of energy among the neutrino 
flavors. After weak decoupling, the v e s and i7 e s remain in equilibrium with the 
matter for a longer period than their heavy-flavor counterparts, due to the 
larger cross sections for scattering off electrons and because of the charge- 
current reactions 

i> e + n «-> p + e~ 

i7 e + p <-> n + e + . (71) 

Thus the heavy flavor neutrinos decouple from deeper within the star, where 
temperatures are higher. Typical calculations yield 

T Vfl ~ T Vt ~ 8MeV T„ e ~ 3.5MeV T„- - 4.5MeV. (72) 

The difference between the v e and v e temperatures is a result of the neutron 
richness of the matter, which enhances the rate for charge-current reactions of 
the v e s, thereby keeping them coupled to the matter somewhat longer. 

This temperature hierarchy is crucially important to nucleosynthesis and 
also to possible neutrino oscillation scenarios. The three- flavor MSW level- 
crossing diagram is shown in Fig. 10. One very popular scenario attributes 
the solar neutrino problem to <-> v e transmutation; this means that a sec- 
ond crossing with a v T could occur at higher density. It turns out plausible 
seasaw mass patterns suggest a v T mass on the order of a few eV, which would 
be interesting cosmologically. The second crossing would then occur outside 
the neutrino sphere, that is, after the neutrinos have decoupled and have fixed 
spectra with the temperatures given above. Thus a v e <-> v T oscillation would 
produce a distinctive T ~ 8 MeV spectrum of v e s. This has dramatic conse- 
quences for terrestrial detection and for nucleosynthesis in the supernova. 

4.2 The Neutrino Process® 

Core-collapse supernovae are one of the major engines driving galactic chem- 
ical evolution, producing and ejecting the metals that enrich our galaxy. The 
discussion of the previous section described the hydrostatic evolution of a pre- 
supernova star in which large quantities of the most abundant metals (C, O, Ne, 
...) are synthesized and later ejected during the explosion. During the passage 
of the shock wave through the star's mantle, temperature of ~ (1 — 3) • 10 9 K and 
are reached in the silicon, oxygen, and neon shells. This shock wave heating 
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~10 g/cm density vacuum 

Figure 10: Three-flavor neutrino level-crossing diagram. One popular scenario associates the 
solar neutrino problem with v £ oscillations and predicts a cosmologically interested 

massive v T with v e <-> u T oscillations near the supernova ncutrinosphere. 

induces (7, a) <-» (a, 7) and related reactions that generate a mass flow toward 
highly bound nuclei, resulting in the synthesis of iron peak elements as well as 
less abundant odd-A species. Rapid neutron-induced reactions are thought to 
take place in the high-entropy atmosphere just above the mass cut, producing 
about half of the heavy elements above A ~ 80. This is the subject of the Sec. 
4.3. Finally, the ^-process described below is responsible for the synthesis of 
rare species such as n B and 19 F. This process involves the response of nuclei 
at momentum transfers where the allowed approximation is no longer valid. 
Thus we will use the ^-process in this section to illustrate some of the relevant 
nuclear physics. 

One of the problems - still controversial - that may be connected with the 
neutrino process is the origin of the light elements Be, B and Li, elements which 
are not produced in sufficient amounts in the big bang or in any of the stellar 
mechanisms we have discussed. The traditional explanation has been cosmic 
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ray spallation interactions with C, O, and N in the interstellar medium. In this 
picture, cosmic ray protons collide with C at relatively high energy, knocking 
the nucleus apart. So in the debris one can find nuclei like 10 B, n B, and 7 Li. 

But there are some problems with this picture. First of all, this is an 
example of a secondary mechanism: the interstellar medium must be enriched 
in the C, O, and N to provide the targets for these reactions. Thus cosmic 
ray spallation must become more effective as the galaxy ages. The abundance 
of boron, for example, would tend to grow quadratically with metalicity, since 
the rate of production goes linearly with metalicity. But observations, espe- 
cially recent measurements with the HST, find a linear growthEO in the boron 
abundance. 

A second problem is that the spectrum of cosmic ray protons peaks near 
1 GeV, leading to roughly comparable production of the two isotopes 10 B 
and 11 B. That is, while it takes more energy to knock two nucleons out of 
carbon than one, this difference is not significant compared to typical cosmic 
ray energies. More careful studies lead to the expectation that the abundance 
ratio of n B to 10 B might be ~ 2. In nature, it is greater than 4. 

Fans of cosmic ray spallation have offered solutions to these problems, 
e.g., similar reactions occurring in the atmospheres of nebulae involving lower 
energy cosmic rays. As this suggestion was originally stimulated by the ob- 
servation of nuclear 7 rays from Orion, now retracted, some of the motivation 
for this scenario has evaporated. Here I focus on an alternative explanation, 
synthesis via neutrino spallation. 

Previously we described the allowed Gamow- Teller (spin-flip) and Fermi 
weak interaction operators. These are the appropriate operators when one 
probes the nucleus at a wavelength - that is, at a size scale - where the nucleus 
responds like an elementary particle. We can then characterize its response by 
its macroscopic quantum numbers, the spin and charge. On the other hand, 
the nucleus is a composite object and, therefore, if it is probed at shorter length 
scales, all kinds of interesting radial excitations will result, analogous to the 
vibrations of a drumhead. For a reaction like neutrino scattering off a nucleus, 
the full operator involves the additional factor 



ik-r 



1 + ik ■ f 



(73) 



where the expression on the right is valid if the magnitude of k is not too large. 
Thus the full charge operator includes a "first forbidden" term 



A 




(74) 
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and similarly for the spin operator 



A 



^2[n <g> a(i)]j =0 ,i,2T 3 (i). 



(75) 



These operators generate collective radial excitations, leading to the so-called 
"giant resonance" excitations in nuclei. The giant resonances are typically at 
an excitation energy of 20-25 MeV in light nuclei. One important property is 
that these operators satisfy a sum rule (Thomas-Reiche-Kuhn) of the form 



where the sum extends over a complete set of final nuclear states. These 
first-forbidden operators tend to dominate the cross sections for scattering the 
high energy supernova neutrinos (v^s and z^s), with E v ~ 25 MeV, off light 
nuclei. From the sum rule above, it follows that nuclear cross sections per 
target nucleon are roughly constant. 

The El giant dipole mode described above is depicted qualitatively in Fig. 
11a. This description, which corresponds to an early model of the giant res- 
onance response by Goldhaber and Teller, involves the harmonic oscillation 
of the proton and neutron fluids against one another. The restoring force for 
small displacements would be linear in the displacement and dependent on the 
nuclear symmetry energy. There is a natural extension of this model to weak 
interactions, where axial excitations occur. For example, one can envision a 
mode similar to that of Fig. 11a where the spin-up neutrons and spin-down pro- 
tons oscillate against spin-down neutrons and spin-up protons, the spin-isospin 
mode of Fig. lib. This mode is one that arises in a simple SU(4) extension of 
the Goldhaber- Teller model, derived by assuming that the nuclear force is spin 
and isospin independent, at the same excitation energy as the El mode. In 
full, the Goldhaber- Teller model predicts a degenerate 15-dimcnsional super- 
multiplet of giant resonances, each obeying sum rules analogous to the TRK 
sum rule. While more sophisticated descriptions of the giant resonance region 
are available, of course, this crude picture is qualitatively accurate. 

This nuclear physics is important to the ^-process. The simplest example 
of ^-process nucleosynthesis involves the Ne shell in a supernova. Because of 
the first-forbidden contributions, the cross section for inelastic neutrino scat- 
tering to the giant resonances in Ne is <~ 3 • 10~ 41 cm 2 /flavor for the more 
energetic heavy-flavor neutrinos. This reaction 




(76) 



v + A -» v' + A* 



(77) 
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Figure 11: Schematic illustration of a) the El giant dipole mode familiar from electromag- 
netic interactions and b) a spin-isospin giant dipole mode associated with the first-forbidden 
weak axial response. 



transfers an energy typical of giant resonances, ~ 20 MeV. A supernova releases 
about 3 x 10 53 ergs in neutrinos, which converts to about 4 x 10 57 heavy flavor 
neutrinos. The Ne shell in a 20 M Q star has at a radius ~ 20,000 km. Thus 
the neutrino fluence through the Ne shell is 

4 • 10 57 

<t> To ~ 10 38 /cm 2 . (78) 

V 4tt(20, 000km) 2 1 1 ; 

Thus folding the fluence and cross section, one concludes that approximately 
1 /300th of the Ne nuclei interact. 

This is quite interesting since the astrophysical origin of 19 F had not been 
understood. The only stable isotope of fluorine, 19 F has an abundance 

19F 1 

(79) 



20 Nc 3100' 

This leads to the conclusion that the fluorine found in toothpaste was created 
by neutral current neutrino reactions deep inside some ancient supernova. 

The calculation of the final 19 F/ 20 Ne ratio is more complicated than the 
simple 1/300 ratio given above: 

• When Ne is excited by <~ 20 MeV through inelastic neutrino scattering, it 
breaks up in two ways 

20 Ne(is, z/) 20 Ne* _»19 Ne + n _>i9 p + c + + „ e + n 
20 Nc(^,^) 2 °Ne* ^ 19 F + p (80) 
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with the first reaction occurring half as frequently as the second. As both 
channels lead to 19 F, we have correctly estimated the instantaneous abundance 
ratio in the Ne shell of 

i9 F 1 

20Ne"~ 300' ^ 

• We must also address the issue of whether the produced 19 F survives. In the 
first 10 -8 sec the coproduced neutrons in the first reaction react via 

15 0(n,p) 15 N 19 Nc(n,a) 16 20 Nc(n, 7 ) 21 Ne 19 Ne(n,p) 19 F (82) 

with the result that about 70% of the 19 F produced via spallation of neutrons 
is then immediate destroyed, primarily by the (n,a) reaction above. In the 
next 10 -6 sec the coproduced protons are also processed 

15 N(p,a) 12 C 19 F(p,a) 16 23 Na(p, a) 20 Ne (83) 

with the latter two reactions competing as the primary proton poisons. This 
makes an important prediction: stars with high Na abundances should make 
more F, as the 23 Na acts as a proton poison to preserve the produced F. 

• Finally, there is one other destruction mechanism, the heating associated 
with the passage of the shock wave. It turns out the the F produced prior to 
shock wave passage can survive if it is in the outside half of the Ne shell. The 
reaction 

19 F( 7 ,a) 15 N (84) 

destroys F for peak explosion temperatures exceeding 1.7-10 9 K. Such a temper- 
ature is produced at the inner edge of the Ne shell by the shock wave heating, 
but not at the outer edge. 

If all of this physics in handled is a careful network code that includes 
the shock wave heating and F production both before and after shock wave 
passage, the following are the results: 

[ 19 F/ 20 Ne]/[ 19 F/ 20 Ne] Q T hcavy „(MeV) 



0.14 4 

0.6 6 

1.2 8 

1.1 10 

1.1 12 



where the abundance ratio in the first column has been normalized to the solar 
value. One sees that the attribution of F to the neutrino process argues that 
the heavy flavor v temperature must be greater than 6 MeV, a result theory 
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favors. One also sees that F cannot be overproduced by this mechanism: 
although the instantaneous production of F continues to grow rapidly with 
the neutrino temperature, too much F results in its destruction through the 
(p, a) reaction, given a solar abundance of the competing proton poison 23 Na. 
Indeed, this illustrates an odd quirk: although in most cases the neutrino 
process is a primary mechanism, one needs 23 Na present to produce significant 
F. Thus in this case the neutrino process is a secondary mechanism. 

While there are other significant neutrino process products ( 7 Li, 138 La, 
180 Ta, 15 N ...), the most important product is n B, produced by spallation off 
carbon. A calculation by Timmes et al. [18] found that the combination of the 
neutrino process, cosmic ray spallation and big-bang nucleosythesis together 
can explain the evolution of the light elements. The neutrino process, which 
produces a great deal of 11 B but relatively little 10 B, combines with the cos- 
mic ray spallation mechanism to yield the observed isotope ratio. Again, one 
prediction of this picture is that early stars should be 11 B rich, as the neutrino 
process is primary and operates early in our galaxy's history; the cosmic ray 
production of 10 B is more recent. There is hope that HST studies will soon be 
able to descriminate between 10 B and n B: as yet this has not been done. 

4-3 The r-process 

Beyond the iron peak nuclear Coulomb barriers become so high that charged 
particle reactions become ineffective, leaving neutron capture as the mechanism 
responsible for producing the heaviest nuclei. If the neutron abundance is 
modest, this capture occurs in such a way that each newly synthesized nucleus 
has the opportunity to (3 decay, if it is energetically favorable to do so. Thus 
weak equilibrium is maintained within the nucleus, so that synthesis is along 
the path of stable nuclei. This is called the s- or slow-process. However a plot 
of the s-process in the (N,Z) plane reveals that this path misses many stable, 
neutron-rich nuclei that are known to exist in nature. This suggests that 
another mechanism is at work, too. Furthermore, the abundance peaks found 
in nature near masses A <~ 130 and A <~ 190, which mark the closed neutron 
shells where neutron capture rates and (3 decay rates are slower, each split into 
two subpeaks. One set of subpeaks corresponds to the closed-neutron-shell 
numbers N <~ 82 and N ~ 126, and is clearly associated with the s-process. The 
other set is shifted to smaller N, <~ 76 and ~ 116, respectively, and is suggestive 
of a much more explosive neutron capture environment where neutron capture 
can be rapid. 

This second process is the r- or rapid-process, characterized by: 
• The neutron capture is fast compared to f3 decay rates. 
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• The equilibrium maintained within a nucleus is established by (n, 7) <-> (7, n): 
neutron capture fills up the available bound levels in the nucleus until this 
equilibrium sets in. The new Fermi level depends on the temperature and the 
relative n/7 abundance. 

• The nucleosynthesis rate is thus controlled by the (3 decay rate: each (3~ 
capture coverting n — > p opens up a hole in the neutron Fermi sea, allowing 
another neutron to be captured. 

• The nucleosynthesis path is along exotic, neutron-rich nuclei that would be 
highly unstable under normal laboratory conditions. 

• As the nucleosynthesis rate is controlled by the (3 decay, mass will build up 
at nuclei where the (3 decay rates are slow. It follows, if the neutron flux is 
reasonable steady over time so that equilibrated mass flow is reached, that the 
resulting abundances should be inversely proportional to these (3 decay rates. 

Let's first explore the (n, 7) <-> (7, n) equilibrium condition, which requires 
that the rate for (71,7) balances that for (7,71) for an average nucleus. So 
consider the formation cross section 

A + n^{A+l)+ 1 . (85) 

This is an exothermic reaction, as the neutron drops into the nuclear well. Our 
averaged cross section, assuming a resonant reaction (the level density is high 
in heavy nuclei) is 

where E <~ is the resonance energy, and the Ts are the indicated partial and 
total widths. Thus the rate per unit volume is 

where N n and Na are the neutron and nuclear number densities and // the 
reduced mass. This has to be compared to the (7, n) rate. 

The (7, n) reaction requires the photon number density in the gas. This 
is given by the Bose-Einstein distribution 

The high-energy tail of the normalized distribution can thus be written 



l —e 2 e-^ kT de (89) 



iV 7 7T 2 
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where in the last expression we have set h = c = 1. 

Now we need the resonant cross section in the (7, n) direction. For photons 
the wave number is proportional to the energy, so 



7r r 7 r„ 
T(7 ' n) ~ ^ ( e _ Er f + (r/2) 2 ' 



(90) 



As the velocity is c =1, 



^-Ts-S^^ ^wm - (91) 

We evaluate this in the usual way for a sharp resonance, remembering that the 
energy integral over just the denominator above (the sharply varying part) is 
2tt/T 

~ T -^e- E ^_ (92) 



So that the rate becomes 

r {l>n) ~2N A+1 T -^e- E rl kT . (93) 
Equating the (n, 7) and (7, n) rates and taking Na ~ AU-i then yields 



Nn-j^A^-X' 2 e~ E ^ kT (94) 



(?ic) 3 V 2tt 

where the hs and cs have been properly inserted to give the right dimensions. 
Now E r is esssentially the binding energy. So plugging in the conditions N n ~ 
3 x 10 23 /cm 3 and Tg ~ 1, we find that the binding energy is ~ 2.4 MeV. Thus 
neutrons are bound by about 30 times kT, a value that is still small compared 
to a typical binding of 8 MeV for a normal nucleus. (In this calculation I 
calculated the neutron reduced mass assuming a nuclear target with A=150.) 

The above calculation fails to count spin states for the photons and nuclei 
and is thus not quite correct. But it makes the essential point: the r-process 
involves very exotic species largely unstudied in any terrestrial laboratory. It 
is good to bear this in mind, as in the following section we will discuss the 
responses of such nuclei to neutrinos. Such responses thus depend on the 
ability of theory to extrapolate responses from known nuclei to those quite 
unfamiliar. 

The path of the r-process is along neutron-rich nuclei, where the neutron 
Fermi sea is just ~ (2-3) MeV away from the neutron drip line (where no 
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more bound neutron levels exist). After the r- process finishes (the neutron 
exposure ends) the nuclei decay back to the valley of stability by (3 decay. This 
can involve some neutron spallation (/3-delayed neutrons) that shift the mass 
number A to a lower value. But it certainly involves conversion of neutrons 
into protons, and that shifts the r-process peaks at N ~ 82 and 126 to a lower 
N, off course. This effect is clearly seen in the abundance distribution: the 
r-process peaks are shifted to lower N relative to the s-process peaks. This is 
the origin of the second set of "subpeaks" mentioned at the start of the section. 

It is believed that the r-process can proceed to very heavy nuclei (A ~ 
270) where it is finally ended by /3-delayed and n-induced fission, which feeds 
matter back into the process at an A ~ A max /2. Thus there may be important 
cycling effects in the upper half of the r-process distribution. 

What is the site(s) of the r-process? This has been debated many years 
and still remains a controversial subject: 

• The r-process requires exceptionally explosive conditions 

p(n) ~ 10 20 cm" 3 T - 10 9 K t - Is. 

• Both primary and secondary sites proposed. Primary sites are those not 
requiring preexisting metals. Secondary sites are those where the neutron cap- 
ture occurs on preexisting s-process seeds. 

• Suggested primary sites include the the neutronized atmosphere above the 
proto-neutron star in a Type II supernova, neutron-rich jets produced in su- 
pernova explosions or in neutron star mergers, inhomogeneous big bangs, etc. 

• Secondary sites, where p(n) can be lower for successful synthesis, include the 
He and C zones in Type II supernovae, the red giant He flash, etc. 

The balance of evidence favors a primary site, so one requiring no preen- 
richment of heavy s-process metals. Among the evidence: 

1) HST studies of very- metal-poor halo stars: The most important evidence 
are the recent HST measurements of Cowan, Sneden et al. E3 of very metal- 
poor stars ([Fe/H] ~ -1.7 to -3.12) where an r-process distribution very much 
like that of our sun has been seen for Z > 56. Furthermore, in these stars 
the iron content is variable. This suggests that the "time resolution" inher- 
ent in these old stars is short compared to galactic mixing times (otherwise 
Fe would be more constant). The conclusion is that the r-process material 
in these stars is most likely from one or a few local supernovae. The fact 
that the distributions match the solar r-process (at least above charge 56) 
strongly suggests that there is some kind of unique site for the r-process: the 
solar r-process distribution did not come from averaging over many different 
kinds of r-process events. Clearly the fact that these old stars are enriched 
in r-process metals also strongly argues for a primary process: the r-process 
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works quite well in an environment where there are few initial s-process metals. 

2) There are also fairly good theoretical arguments that a primary r-process oc- 
curring in a core-collapse supernova might be viableE3. First, galactic chemical 
evolution studies indicate that the growth of r-process elements in the galaxy 
is consistent with low-mass Type II supernovae in rate and distribution. More 
convincing is the fact that modelers have shown that the conditions needed for 
an r-process (very high neutron densities, temperatures of 1-3 billion degrees) 
might be realized in a supernova. The site is the last material expelled from 
the supernova, the matter just above the mass cut. When this material is 
blown off the star initially, it is a very hot neutron-rich, radiation-dominated 
gas containing neutrons and protons, but an excess of the neutrons. As it 
expands off the star and cools, the material first goes through a frcezeout to 
a particles, a step that essentially locks up all the protons in this way. Then 
the as interact through reactions like 

a + a + a -> 12 C 
a + a + n — > 9 Be 

to start forming heavier nuclei. Note, unlike the big bang, that the density 
is high enough to allow such three-body interactions to bridge the mass gaps 
at A = 5,8. The a capture continues up to heavy nuclei, to A ~ 80, in the 
network calculations. The result is a small number of "seed" nuclei, a large 
number of as, and excess neutrons. These neutrons preferentially capture on 
the heavy seeds to produce an r-process. Of course, what is necessary is to 
have ~ 100 excess neutrons per seed in order to successfully synthesize heavy 
mass nuclei. Some of the modelers find conditions where this almost happens. 

There are some very nice aspects of this site: the amount of matter ejected 
is about 10~ 5 — 10~ 6 solar masses, which is just about what is needed over the 
lifetime of the galaxy to give the integrated r-process metals we see, taking a 
reasonable supernova rate. But there are also a few problems, especially the 
fact that with calculated entropies in the nucleon soup above the proto-neutron 
star, neutron fractions appear to be too low to produce a successful A ~ 190 
peak. There is some interesting recent work invoking neutrino oscillations 
to cure this problem: charge current reactions on free protons and neutrons 
determine the n/p ratio in the gas. Then, for example, an oscillation of the 
type v e — > I'steriie can alter this ratio, as it would turn off the v e s that destroy 
neutrons by charged-current reactions. Unfortunately, a full discussion of such 
possibilities would take us too far afield today. 

The nuclear physics of the r-process tells us that the synthesis occurs when 
the nucleon soup is in the temperature range of (3-1) -10 9 K, which, in the hot 
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bubble r-process described above, corresponds to a freezeout radius of (600- 
100) km and a time ~ 10 seconds after core collapse. The neutrino fluence after 
freezeout (when the temperature has dropped below 10 9 K and the r-process 
stops) is then ~ (0.045-0.015) -10 51 ergs/(100km). Thus, after completion of 
the r-process, the newly synthesized material experiences an intense flux of 
neutrinos. This brings up the question of whether the neutrino flux could have 
any effect on the r-process. 

4-4 Neutrinos and the r-process^ 

Rather than describe the exotic effects of neutrino oscillations on the r-process, 
mentioned briefly above, we will examine standard-model effects that are never- 
theless quite interesting. The nuclear physics of this section neutrino-induced 
neutron spallation reactions - is also relevant to recently proposed supernova 
neutrino observatories such as OMNIS and LAND. In contrast to our first 
discussion of the z/-process in Sec. 4.2, it is apparent that neutrino effects 
could be much larger in the hot bubble r-process: the synthesis occurs much 
closer to the star than our Ne radius of 20,000 km: estimates are 600-1000 
km. The r-process is completed in about 10 seconds (when the temperature 
drops to about one billion degrees), but the neutrino flux is still significant as 
the r-process freezes out. The net result is that the "post-processing" neutrino 
fluence - the fluence that can alter the nuclear distribution after the r-process 
is completed - is about 100 times larger than that responsible for fluorine 
production in the Ne zone. Recalling that 1/300 of the nuclei in the Ne zone 
interacted with neutrinos, and remembering that the relevant neutrino- nucleus 
cross sections scale as A, one quickly sees that the probability of a r-process 
nucleus interacting with the neutrino flux is approximately unity. 

Because the hydrodynamic conditions of the r-process are highly uncer- 
tain, one way to attack this problem is to work backward in time. We know 
the final r-process distribution (what nature gives us) and we can calculate 
neutrino-nucleus interactions relatively well. Thus from the observed r-process 
distribution (including neutrino postprocessing) we can work backward to find 
out what the r-process distribution looked like at the point of freezeout. In 
Figs. 12 and 13, the "real" r-process distribution - that produced at freeze- 
out - is given by the dashed lines, while the solid lines show the effects of the 
neutrino postprocessing for a particular choice of fluence. The nuclear physics 
input into these calculations is precisely that previously described: GT and 
first-forbidden cross sections, with the responses centered at excitation ener- 
gies consistent with those found in ordinary, stable nuclei, taking into account 
the observed dependence on \N — Z\. 



52 





110 115 120 



125 130 135 140 
A 



Figure 12: Comparison of the r-process distribution that would result from the freezeout 
abundances near the A ~ 130 mass peak (dashed line) to that where the effects of neutrino 
postprocessing have been include (solid line). The fluence has been fixed by assuming that 
the A = 124-126 abundances are entirely due to the ^-process. 

One important aspect of the figures is that the mass shift is significant. 
This has to do with the fact that a 20 MeV excitation of a neutron-rich nucleus 
allows multiple neutrons ( <~ 5) to be emitted. (Remember we found that the 
binding energy of the last neutron in an r-proccss neutron- rich nuclei was about 
2-3 MeV under typical r-process conditions.) The second thing to notice is that 
the relative contribution of the neutrino process is particularly important in 
the "valleys" beneath the mass peaks: the reason is that the parents on the 
mass peak are abundant, and the valley daughters rare. In fact, it follows 
from this that the neutrino process effects can be dominant for precisely seven 
isotopes (Te, Re, etc.) lying in these valleys. Furthermore if an appropriate 
neutrino fluence is picked, these isotope abundances are produced perfectly 
(given the abundance errors). The fluences are 

N = 82 peak 0.031 • 10 51 ergs/(100km) 2 /flavor 



•53 



N = 126 peak 0.015 • 10 51 crgs/(100km) 2 /flavor, 

values in fine agreement with those that would be found in a hot bubble r- 
process. So this is circumstantial but significant evidence that the material 
near the mass cut of a Type II supernova is the site of the r-process: there is 
a neutrino fingerprint. 

.6 



.5 



.4 

.2 



.1 





170 175 180 185 190 195 200 
A 

Figure 13: As in Fig. 12, but for the A ~ 195 mass peak. The A = 183-187 abundances are 
entirely attributed to the (/-process. 



4-5 Neutrino Oscillations and the r-process 

For the usual seesaw pattern of neutrino masses and a cosmological interesting 
v T (i.e., a heavy neutrino with a mass in the neighborhood of 10 cV), the full 
MSW pattern is shown in Fig. 10. If the v e — crossing is responsible for 
the solar neutrino problem, a second crossing, v e — v Tl is expected at a density 
large compared to that of the solar core, but small compared to the location 
of the supernova neutrinosphere (~ 10 12 g/cm 3 ). For a very large range of 
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mixing angles, this crossing is adiabatic and thus leads to v e <-» v T conversion. 
These spectra thus change identities, leading to an anomalously hot v e flux 
from a Type II supernova. 

As the i^-nucleon cross section is proportional to E 2 , the reaction v e + n 
— ► e~ + p is enhanced, while v e + p — * c + + n is unchanged. For a rather 
extensive range of v e <-> v T mixing angles and Sm 2 , this crossing then destroys 
the r-process: the hotter v e s drive the matter proton richllJ. Thus, if one 
accepts this location as the site of the r-process, very strong constraints on 
cosmologically interesting v T s are obtained. These limits are truly remarkable 
for their sensitivity to small mixing angles, extended to sin 2 29 ~ 10~ 5 for 
neutrino mass differences above a few eV 2 . 

I thank A. S. Brun and R. E. Shrock for helpful comments, and Paul 
Langacker for his able organization of the 1998 TASI summer school. This 
work was supported in part by the US Department of Energy. 



References 

1. This section taken from W. C. Haxton, Ann. Rev. Astron. Astrophys. 
33, 459 (1995). 

2. R. Davis, Jr., D. S. Harmer, and K. C. Hoffman, Phys. Rev. Lett. 20, 
1205 (1966). 

3. J. N. Abdurashitov et al., Phys. Lett. B 328, 234 (1994). 

4. P. Anselmann et al., Phys. Lett. B 285, 376 (1992). 

5. Y. Suzuki, Nucl. Phys. B 38, 54 (1995). 

6. Y. Suzuki, talk presented at Neutrino '98 (Takayama, Japan, June, 1998). 

7. T. Kajita, talk presented at Neutrino '98 (Takayama, Japan, June, 1998). 

8. J. N. Bahcall, S. Basu, and M. H. Pinsonneault, Phys. Lett. B 433, 1 
(1998). 

9. A. S. Brun, S. Turck-Chieze, and P. Morel, Ap. J. 506, 913 (1998) and 
private communication; S. Turck-Chieze and I. Lopez, Ap. J. 408, 347 
(1993). 



10. This section taken from W. C. Haxton, |nucl-th/9901037 . 

11. S. M. Austin, N. Anantaraman, and W. G. Love, Phys. Rev. Lett. 73, 
30 (1994); J. W. Watson et al., Phys. Rev. Lett. 55, 1369 (1985). 

12. K. Lande, talk presented at Neutrino '98 (Takayama, Japan, June, 1998). 

13. J. N. Bahcall, Neutrino Astrophysics, (Cambridge University, Cam- 
bridge, 1989). 

14. C. W. Johnson, E. Kolbe, S. E. Koonin, and K. Langanke, Ap. J. 392, 
320 (1992). 

15. B. W. Filippone, A. J. Elwyn, C. N. Davids, and D. D. Koetke, Phys. 



55 



Rev. Lett. 50, 412 (1983). 

16. J. N. Bahcall and R. Ulrich, Rev. Mod. Phys. 60, 297 (1988). 

17. J. N. Bahcall and W. C. Haxton, Phys. Rev. D 40, 931 (1989). 

18. V. Castellani, S. Degl'Innocenti, G. Fiorentini, M. Lissia, and B. Ricci, 
Phys. Rev. D 50, 4749 (1994). 

19. J. N. Bahcall and R. Davis, Jr., in Essays in Nuclear Astrophysics, ed. 

C. A. Barnes, D. D. Clayton, and D. Schramm (Cambridge Univ. Press, 
Cambridge) p. 243. 

20. F. W. W. Dilke and D. O. Cough, N 240, 262 (1972). 

21. W. J. Merryfield, in Solar Modeling, ed. A. B. Balantekin and J. N. 
Bahcall (World Scientific, Singapore, 1995). 

22. A. Cumming and W. C. Haxton, Phys. Rev. Lett. 77, 4286 (1996). 

23. N. Hata and P. Langacker, Phys. Rev. D 56, 6107 (1997). 

24. K. M. Heeger and R. G. H. Robertson, Phys. Rev. Lett. 77, 3720 (1996). 

25. For a proper (wave packet) derivation see M. Nauenberg, submitted to 
Phys. Lett. B. 

26. S. P. Mikheyev and A. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985); L. 
Wolfenstein, Phys. Rev. D 17, 2369 (1979). 

27. H. Bethe, Phys. Rev. Lett. 56, 1305 (1986). 

28. W. C. Haxton, Phys. Rev. Lett. 57, 1271 (1986). 

29. S. J. Parke, Phys. Rev. Lett. 57, 1275 (1986). 

30. C. S. Lim and W. J. Marciano, Phys. Rev. D 37, 1368 (1988). 

31. E. Kh. Akhmedov, Sov. J. Nucl. Phys. 48, 382 (1988). 

32. G. Raffelt, Phys. Rev. Lett. 64, 2856 (1990). 

33. K. Fujikawa and R. E. Shrock, Phys. Rev. Lett. 43, 963 (1980). 

34. W. C. Haxton and G. J. Stephenson, Jr., Prog. Part. Nucl. Phys. 12, 
409 (1984). 

35. A. Mezzacappa et al., Ap. J. 495, 911 (1998); H.-Th. Janka and E. 
Muller, Astron. Astrophys. 306, 167 (1996); A. Burrows, S. Hayes, and 
B. A. Fryxell, Ap. J. 450, 830 (1995). 

36. S. E. Woosley and W. C. Haxton, Nature 334, 45 (1988); S. E. Woosley, 

D. H. Hartmann, R. D. Hoffman, and W. C. Haxton, Ap. J. 356, 272 
(1990). 

37. F. X. Timmes, S. E. Woosley, and T. A. Weaver, Ap. J. Suppl. 98, 617 
(1995). 

38. J. J. Cowan et al., astro-ph/9808272 (to appear in Ap. J). 

39. S. E. Woosley J. R. Wilson, G. J. Mathews, R. D. Hoffman, and B. S. 
Meyer, Ap. J. 433, 229 (1994). 

40. W. C. Haxton, K. Langanke, Y.-Z. Qian, and P. Vogel, Phys. Rev. Lett. 
78, 2694 (1997) and Phys. Rev. C 55, 1532 (1997). 



56 



41. Y.-Z. Qian ct al., Phys. Rev. Lett. 71, 1965 (1993). 



57 



